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ABSTRACT 

Using a model for the self-regulated growth of supermassive black holes in mergers involving gas-rich galax- 
ies, we study the relationship between quasars and the population of merging galaxies and thereby predict the 
merger-driven star formation rate density of the Universe. In our picture, mergers drive gas inflows, fueling 
nuclear starbursts and "buried" quasars until feedback disperses the gas, allowing the quasar to be briefly visi- 
ble as a bright optical source. As black hole accretion declines, the quasar dies and the stellar remnant relaxes 
passively with properties and correlations typical of red, elliptical galaxies. By simulating the evolution of 
such events, we demonstrate that the observed statistics of merger rates/fractions, luminosity functions, mass 
functions, star formation rate distributions, quasar luminosity functions, quasar host galaxy luminosity func- 
tions, and elliptical/red galaxy luminosity and mass functions are self-consistent. We use our simulations to 
de-convolve both the quasar and merging galaxy luminosity functions to determine the birthrate of black holes 
of a given final mass and merger rates as a function of the total stellar mass and the mass of new stars formed 
during a merger. From this, we predict the merging galaxy luminosity function in various observed wavebands 
(e.g. UV, optical, and near-IR), color-magnitude relations, mass functions, absolute and specific star formation 
rate distributions and star formation rate density, and quasar host galaxy luminosity function, as a function of 
redshift from z = 0-6. We invert this relationship to predict e.g. quasar luminosity functions from observed 
merger luminosity functions or star formation rate distributions. Our results show good agreement with obser- 
vations, but idealized models of quasar lightcurves give inaccurate estimates and are ruled out by comparison of 
merging galaxy and quasar observations at > 99.9% confidence, provided that quasars are triggered in mergers. 
Using only observations of quasars, we estimate the contribution of mergers to the star formation rate density 
of the Universe to high redshifts, z ~ 4, and constrain the evolution in the characteristic initial gas fractions of 
quasar and spheroid-producing mergers. 

Subject headings: quasars: general — galaxies: active — galaxies: evolution — cosmology: theory 



1. INTRODUCTION 

It is believed that our Universe evolved hierarchically, with 
small objects merging to form progressively larger systems. 
Observations of interacting galaxies support this view and led 
Toomre & Toomre (1972) and Toomre (1977) to propose the 
"merger hypothesis," according to which ellipticals originate 
when spiral galaxies collide (for reviews, see e.g. Barnes & 
Hernquist 1992; Barnes 1998; Schweizer 1998). It has also 
been established that supermassive black holes are ubiquitous 
in the centers of galaxies (e.g. Kormendy & Richstone 1995; 
Richstone et al. 1998; Kormendy & Gebhardt 2001) and that 
their masses correlate with either the masses (Magorrian et 
al. 1998) or the velocity dispersions (i.e. the Mbh-c rela- 
tion: Ferrarese & Merritt 2000; Gebhai'dt et al. 2000) of their 
host spheroids. These relations argue that supermassive black 
holes and galaxies are linked in an evolutionary manner, as 
implied by simulations showing that self-regulated black hole 
growth in a galaxy merger has a significant impact on the 
structure of the remnant (Di Matteo et al. 2005). 

Observations also indicate that galaxy mergers produce 
starbursts and fuel black hole growth, and that both play roles 
in structuring ellipticals. Infrared luminous galaxies appear 
to be powered at least in part by nuclear starbursts (e.g. Soifer 
et al. 1984a,b; Sanders et al. 1986, 1988a,b; for a review, 
see e.g. Soifer et al. 1987), and the most intense of these, 
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ultraluminous infrared galaxies (ULIRGs), are always associ- 
ated with mergers (e.g. Allen et al. 1985; Joseph & Wright 
1985; Armus et al. 1987; Kleinmann et al. 1988; Melnick & 
Mirabel 1990; for reviews, see Sanders & Mirabel 1996 and 
Jogee 2005). Radio observations reveal large quantities of 
dense gas in the centers of ULIRGs (e.g. Scoville et al. 1986; 
Sargent et al. 1987, 1989), providing material to feed black 
hole growth and to boost the concentration and central phase 
space density to match those of ellipticals (see, e.g. Carlberg 
1986; Gunn 1987; Lake 1989; Hernquist 1992, 1993a; Hern- 
quist et al. 1993; Robertson et al. 2006a). Some ULIRGs have 
"warm" IR spectral energy distributions (SEDs), suggesting 
that they harbor buried quasars (e.g. Sanders et al. 1988c), 
as implied by X-ray observations of growing black holes in 
NGC6240 (Komossa et al. 2003) and other ULIRGs, which 
are heavily obscured at optical and infrared wavelengths (e.g. 
Gerssen et al. 2004; Max et al. 2005; Alexander et al. 2005a,b; 
Borys et al. 2005). 

These lines of evidence, and the overlap between the bolo- 
metric luminosities of ULIRGs and quasars, suggest that 
quasars originate from an IR luminous phase of galaxy evolu- 
tion caused by mergers (Sanders et al. 1988a). This proposal 
is supported by observations of nearby hosts of bright quasars 
which reveal features characteristic of mergers (e.g. Stockton 
1978; Heckman et al. 1984; Stockton & MacKenty 1987; 
Stockton & Ridgway 1991; Hutchings & Neff 1992; Bah- 
cafl et al. 1994, 1995, 1997; Canalizo & Stockton 2001) and 
starbursts (e.g., Brotherton et al. 1999. 2002 UCanaUzo et alJ 
■2000; Yip et al. 2004; Vanden Berk et al. 200%. 

Observations of individual merging systems and merger 
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remnants (e.^.. [Lake & Dresslej Il986t iDovon et all |1994|: 
lOliva et aljTi995t LShier & Fischer 1998t LTames etall [19991) 
have shown that their kinematic and photometric properties, 
including velocity dispersions, concentrations, stellar masses, 
light profiles, and phase space densities, are consistent with 
their eventual evolution into typical ~ elliptical galaxies. 
Furthermore, t he correlations obeyed by these mergers and 
remnants (e.g.. iGenzel et al.ll200U iRothberg & Josephll2004 
I2OO6) are similar to e.g. the observed fundamental plane and 
Kormendy relations for relaxed ellipticals, and consistent with 
evolution onto these relations as their stellar populations age. 

While these studies suggest that nearby gas-rich mergers 
pass through starburst and quasar phases and evolve into el- 
lipticals, the implications of these findings at higher redshifts 
are unclear and little is known about the statistics of galaxy 
mergers in this context. In particular, it is difficult to measure 
the detailed properties of d i stant quasar hosts ( e.g., Schweizer 
fT98l lDunlopetalJl2003t iFlovd et alJ l2004l) . and it is not 
clear if there are enough mergers to account for the present 
abundance of spheroids or whether all gas-rich or elliptical- 
producing mergers undergo an infrared luminous starburst 
phase. The gal axy merger rate and i t s evolution has been 
estimated (e.g., 'Xu & Sulentic' '199U iBurkley et alJ Il994t 
Carlberg et al. 1994; Sulentic & Rabacal 1 19941 iKeel & Wu 
1995; Soares et al. 1995; Yee & Ellingson 1995 ; Patton et aL 
1^97; Zenf&Koo 1989; Toledo etal. 1999; Le Fevre et a^ 
2000; Patton et al. 2000, 2002; ConseHce et al. 2003; Xu et al. 
2004; Lin etal. 2004; Bundv et al. 2004; Conselice et aL 
2005t Bell et al. 2006b), but the observational uncertainty 
is large. Moreover, the few observational studies of the 
distribution of quasar host luminosities (Bahcall et al. 1993 
McLure et al. 1999; Falomo et al. 2001; Hamilton et al. 20o3 
Jahnke & Wisot zki 2003; Dunlop et aL .2003; Flovd et all 
2004; Vanden B erk et alJ2006h have not yet determined if this 
distribution is expected (as a subset of the merging galaxy 
population), or whether it requires modes of quasar fuel- 
ing aside from mergers of systems having reservoirs of cold 
gas. Some observations (IStra ughnet aL.2006 ) and theoretical 
analyses (Hopkins et al. 2005^ suggest that the stages most 
likely to be identified as mergers can be quite distinct from the 
merger-driven quasar phase, further complicating such com- 
parisons. 

Indirect evidence that mergers play a role in trans- 
forming galaxies also comes from studies of the stellar 
mass budget contained in morphologically or color-selected 
sub -populations as a function of redshift. For example, 
iBorch et alJ (1200 6) find that the stellar mass in red galax- 
ies (known to be spheroid dominated) grows with time 
from z ^ 1 to the present. Furthermore, measurements 
of galaxy stellar mass functions separated by color or by 
morpholog y (Bundv et al. 2005a b; Franceschini et al. 2006; 
iPannellaet al. 2006) show that the "transition mass," above 
which red galaxies dominate the mass function, increases 
with redshift, following a general trend of "cosmic downsiz- 
ing." As low-mass, red galaxies build up, the mass of the 
largest star-forming galaxies, although less well-constrained, 
decreases correspondingly, further supporting the view that 
star formation is "quenched" in these systems and that they 
move onto the red sequence (Bundv et al. 2005b). By compar- 
ing estimates of the transition or quenching mass to the char- 
acteristic masses of quasars (corresponding to the break in the 
observed quasar luminosity function) and merging galaxies, 
iHopkins et a l. (2006d) have shown that these trace the same 
mass and evolve together over the range < z < 3, imply- 



ing that these objects are linked. This connection motivates 
a detailed analysis to understand the relationship between the 
merger and quasar populations. 

The contribution of mergers to the star formation density of 
the Universe is not well known, especially at high redshifts 
(e.g.. iBrinchmann et al. 1998; Menanteau et al. 2001, 200§, 
iBelletalJ 12001 7" Traditional observational measurements 
which require morphological identification of star-forming 
galaxies are difficult to extend to these redshifts and are likely 
incomplete dWolf et a T 2005). making particularly valuable 
any independent constraints that do not require direct morpho- 
logical information, such as those from quasar populations. 
Even fewer estimates exist for the distribution of star for- 
mation rates in merging systems (e.g.. Bell et al. 2005), and 
whether or not the measurements are consistent with gas-rich 
mergers passing through a starburst phase and building up the 
stellar mass observed in the elliptical galaxy population. 

The significant uncertainties in the relation of ellipti- 
cal galaxy, merger remnant, LIRG/ULIRG, merging (pecu- 
liar/interacting) galaxy, starburst, and quasar populations are 
related to the same underlying issue: Are the luminosity and 
mass distributions of quasars, quasar host galaxies, merging 
galaxies, starbursts, and ellipticals self-consistent? If the ori- 
gin of red/elliptical galaxies primarily involves gas-rich merg- 
ers, with an intermediate starburst phase which is terminated 
by feedback from black hole growth dispersing gas, then these 
populations must track each other at all redshifts. 

A long-standing obstacle to relating these populations has 
been the lack of theoretical models which incorporate the rel- 
evant physical processes. Large-scale cosmological simula- 
tions (e.g., Cen et al. 1994; Zhang et al. 1995; Hernquist et aH 
1996; Katzetal. 1996a b; Navarro et al. 1996; Croft eta!] 
1998 ,1999 .2002: Dave et aL .1999: McD onald et all 
200n[ 120061 iH^Tetal J 1200 1| iViel etalJ 120031 1 2004) and 
semi-analytic models (e.g., Kauffmann & Haehnelt 2003 
Cole et al. 2000; Somerville et al. 2001; Volonteri et al. 20(31 
Wvithe & Loeb 2003; Granato et al. 2004; Some rville et all 

■2004a; Baugh et al. .2005: Croton et al 2006 ) have provided 

a framework for understanding the intergalactic medium 
and large-scale structure formation, but either do not resolve 
or do not follow the detailed effects of star formation, 
supernova feedback, black hole accretion and feedback, and 
dust obscuration, which are critical for inferring the relation 
between e.g. quasar and starburst populations and merger 
rates. Instead, these processes are described by idealized, 
tunable parameterizations rather than being modeled dynam- 
ically. Ideally, the prescriptions should be determined from 
physically motivated models, which have been tested against 
observations. 

Recently, progress has been made in this direction by using 
high-resolution simulations of individual galaxies undergoing 
mergers that incorporate the effects of star formation and pres- 
surization of multi-phase interstellar gas by supernova feed- 
back (Springel & Hernauist 2003a) and feedback from black 
hole growth (Springel et al. 20 05h) . In this manner, Di Mat- 
teo et al. (2005) and Springel et al. (2005a,b) have shown 
that the gas inflows produced by gravitational torques during 
a merger (e.g. Barnes & Hernquist 1991, 1996) both trig- 
ger starbursts (e.g. Mihos & Hernquist 1994a, 1996) and fuel 
rapid black hole growth. During most of this phase, the black 
hole is heavily obscured (Hopkins et al. 2005a), but becomes 
briefly visible as a bright optical quasar when feedback expels 
the surrounding gas. As the gas is heated and dispersed, the 
accretion rate declines, leaving a dead quasar in an ordinary 
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galaxy. The self-regulated nature of this process explains ob- 
served correlatio ns between black hole m ass and properties of 
normal galaxies SDi Matte o et al."200^. as well as the color 
distribution of ellipticals dSprinael et al. 2005a). 

Previously (Hopkins et al. 2005a-d, 2006a-c; Lidz et al. 
2006), we have shown that this time evolution yields a self- 
consistent model for many observed properties of quasars, 
with no dependence on a particular prior cosmological dis- 
tribution. In particular, the simulations imply a new descrip- 
tion of quasar lifetimes where the lifetime of a given source 
depends on both the instantaneous and peak luminosities, so 
that for a given peak luminosity the lifetime increases at lower 
luminosities - i.e. quasars spend most of their lives in phases 
fainter than their peak luminosities. This leads to a differ- 
ent interpretation of the quasar luminosity function than is 
implied by idealized models of quasar lightcurves (in which 
quasars grow/decay in a step function or pure exponential 
manner); namely that the steep, bright end of the quasar lu- 
minosity function traces quasars accreting at high Edding- 
ton ratios near their peak luminosities, but the shallow, faint 
end is dominated by quasars with high peak luminosities but 
which are seen in less luminous states. The distribution of 
quasar birthrates as a function of their peak luminosities (final 
black hole masses), which is directly related to e.g. the gas- 
rich galaxy merger rate as a function of final galaxy mass, is 
peaked at a luminosity /mass corresponding to the break in the 
observed quasar luminosity function. Objects near the peak 
in this distribution dominate the observed faint-end quasar lu- 
min osity function in their f ainter, below peak phases. 

In iHopkins et al.l J2006 c). we used our modeling to pre- 
dict the distribution of remnant elliptical galaxies formed in 
mergers. Because the spheroid stellar mass or velocity dis- 
persion correlates with the final black hole mass in our sim- 
ulations, we can use our model of quasar lifetimes and our 
determination of the birthrate of quasars as a function of their 
peak luminosity (final black hole mass) to infer the rate at 
which spheroids with given properties are formed in merg- 
ers as a function of e.g. mass, velocity dispersion, size, and 
redshift. In this manner, we reproduce spheroid mass, lumi- 
nosity, and velocity dispersion functions, age distributions, 
mass-radius and luminosity-size relations, mass-to-light ra- 
tios, color-magnitude relations, and distributions of young 
(blue) spheroids at redshifts z = 0-6. The co-evolution of 
star formation, black hole growth, and obscuration is a key 
element in this analysis, and corresponding predictions made 
either neglecting the role of black hole feedback in terminat- 
ing star formation or by modeling the quasar lightcurve in an 
idealized manner are inaccurate. 

In principle, the merger hypothesis provides a framework 
for describing the co-formation of quasars and spheroids that 
reproduces a wide range of observations. To investigate 
this further, here we examine the statistics and properties of 
merging systems. In particular, we use our earlier model- 
ing to study the implied relation between the quasar lumi- 
nosity function (QLF), merging galaxy luminosity function 
(MGLF), quasar host galaxy luminosity function (HGLF), 
distribution of star formation rates in mergers (star formation 
rate function; SFRF), and the cumulative merger rates and star 
formation rate density triggered by mergers. 

In §|2j we describe our simulations, and in §|3lwe use them 
to determi ne t he evoluti on o f near-IR and optical/UV lumi- 
nosities (§ 13. U . colors (§ 13. 2t . and star formation rates (§ 13. 3> 
during mergers, as a function of host galaxy properties. In §|3 
we describe our methodology for mapping between various 



distributions (§ 14. 1> . and disc uss t he impact and importance of 
possible systematic effects (§ 14. 2> . In §|5lwe use our modeling 
to predict the quasar luminosity function from the observed 
merging galaxy function (at the same redshift), and vice versa 
(§ 15. 1> . in both /T-band and the optical/UV. We use this to 
predict the underlying birthrate of spheroids/quasars with a 
given final mass/peak luminosity, and compare this with the 
observed merger mass functions (§ 15. 2t . Likewise, we use this 
to predict the distribution of star formation rates in mergers, 
independently from the observed quasar and merger luminos- 
ity functions, and compare with the observed star formation 
rate functions (§ 15. 3> . In § |6l we use the observed evolu- 
tion in the quasar luminosity function to predict the evolution 
of the merging galaxy luminosity function and star formation 
rate function, as well as integrated quantities such as the lu- 
minosity, mass, and star formation rate densities in mergers, 
at redshifts z = 0-6. In §0we consider the relation between 
the MGLF and HGLF, predicting the quasar host galaxy lumi- 
nosity function, the joint distribution of observed quasar and 
quasar host galaxy luminosities, and their co-evolution. Fi- 
nally, in §|8j we summarize our results and discuss their impli- 
cations for the merger hypothesis, various theoretical models, 
and future observations. 

Our approach is described in detail in § I2I4I but readers 
interested in the main results may wish to skip to § |5] or 
§ where we compare and test the relation between vari- 
ous stages and measures of merger-driven activity, and § |6l 
which uses these relations to predict evolution of these mea- 
sures with redshift. 

Throughout, we adopt a JIm = 0.3, fl\ = 0.7, Ho = 
70kms"'Mpc"' cosmology. All magnitudes are Vega mag- 
nitudes, unless otherwise stated. 

2. THE SIMULATIONS 

Our merger simulations were performed with the paral- 
lel TreeSPH code GADGET-2 (Snrinael 2005), b ased on a 
fully conservative formulation ( Sprinael & Hernauist 200^ 
of smoothed particle hydrodynamics (SPH), which conserves 
energy and entropy simultaneously even when smoothing 
lengths evolve adaptively (see e.g., Hernquist 1993b, O'Shea 
et al. 2005). Our simulations account for radiative cooling, 
heating by a UV background (as in Katz et al. 1996, Dave et 
al. 1999), and incorporate a sub-resolution model of a mul- 
tiphase interstellar mediu m (ISM) to describe star for mation 
and supernova feedback ( Soringel & Hernauistl2003al) . Feed- 
back from supernovae is captured in this sub-resolution model 
through an effective equation of state for star-forming gas, en- 
abling us to stably evolve disks with arbitrary gas fractions 
(see, e.g. Springel et al. 2005b; Springel & Hernquist 2005; 
Robertson et al. 2004, 2005). This is described by the param- 
eter qE05, which ranges from ^eos = for an isothermal gas 
with effective temperature of 10"* K, to ^eos = 1 for our full 
multiphase model with an effective temperature 10^ K. 

Supermassive black holes are represented by "sink" par- 
ticles that accrete gas at a rate M estimated from the local 
gas density and sound speed using an Eddington-limited pre- 
scription based on Bondi-Hoyle-Lyttleton accretion theory. 
The bolometric luminosity of the black hole is taken to be 
Lboi = £rMc^, where e,- = 0.1 is the radiative efficiency. We 
assume that a small fraction (typically « 5%) of Lboi couples 
dynamically to the surrounding gas, and that this feedback is 
injected into the gas as thermal energy, weighted by the SPH 
smoothing kernel. This fraction is a free parameter, which 
we determine as in iDi Matteo et alJ J2005h by matching the 
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observed Mbh - c relation. For now, we do not resolve the 
small-scale dynamics of the gas in the immediate vicinity of 
the black hole, but assume that the time-averaged accretion 
rate can be estimated from the gas properties on the scale of 
our spatial resolution (roughly « 20 pc, in the best cases). 

The progenitor galaxy models are described in 
ISpringel et al. (2005b), and we review their properties 
here. For each simulation, we generate two stable, isolated 
di sk galaxies , each with an extended dark matter halo with 
a Hernauist ( 1990) profile, motivated by cosmological 
simulations (e.g. Navarro et al. 1996; Busha et al. 2005), an 
exponential disk of gas and stars, and (optionally) a bulge. 
The galaxies have total masses Mvir = V^jJilOGHo) for z = 0, 
with the baryonic disk having a mass fraction = 0.041, the 
bulge (when present) having m\, = 0.0136, and the rest of the 
mass in dark matter. The dark matter halos are assigned a 
concentration parameter scaled as in Robertson et al. (2006^ 
appropriately for the galaxy mass and redshift following 
[Bullock et al. (2001). The disk scale-length is computed 
based on an assumed spin parameter A = 0.033, chosen to be 
near the mode in the A distribution measured in simulations 
(^itvitska et al. 2002j), and the scale-length of the bulge is set 
to 0.2 times this. 

Typically, each galaxy initially consists of 168000 dark 
matter halo particles, 8000 bulge particles (when present), 
40000 gas and 40000 stellar disk particles, and one BH parti- 
cle. We vary the numerical resolution, with many simulations 
using twice, and a subset up to 128 times, as many particles. 
We choose the initial seed mass of the black hole either in 
accord with the observed Mbh-ct relation or to be sufficiently 
small that its presence will not have an immediate dynamical 
effect, but we have varied the seed mass to identify any sys- 
tematic dependencies. Given the particle numbers employed, 
the dark matter, gas, and star particles are all of roughly equal 
mass, and central cusps in the dark matter and bulge are rea- 
sonably well resolved (see Figure 2 in Springel et al. 2005b). 

Our fitted quasar lifetimes and galaxy scaling relations are 
derived from a series of several hundred simulations of collid- 
ing galaxies, described in detail inlRobertso n et alJ ll2006alhft 
and lHopkins etal] (l2006bV We vary the numerical resolution, 
the orbit of the encounter (disk inclinations, pericenter sepa- 
ration), the masses and structural properties of the merging 
galaxies, initial gas fractions, halo concentrations, the param- 
eters describing star formation and feedback from supernovae 
and black hole growth, and initial black hole masses. 

The progenitor galaxies have virial velocities Vvir = 
80,113,160,226,320, and 500kms-i, and redshifts z = 
0,2, 3, and 6, and our simulations span a range in final black 
hole mass Mbh ^ 10^ - lO'^M©. The extensive range of 
conditions probed provides a large dynamic range, with fi- 
nal spheroid masses spanning ^ 10^- IO'^Mq, covering 
essentially the entire range of the observations we consider at 
all redshifts, and allows us to identify any systematic depen- 
dencies in our models. We consider initial disk gas fractions 
(by mass) of /gas = 0.2, 0.4, 0.8, and 1.0 for several choices 
of virial velocities, redshifts, and ISM equations of state. The 
results described in this paper are based primarily on simula- 
tions of equal-mass mergers; however, by examining a small 
set of simulations of unequal mass mergers, we find that the 
behavior does not change dramatically for mass ratios down to 
about 1:3 or 1:4. This range is appropriate to the observations 
of merging galaxies used in this paper, which are restricted to 
'major' merger events. 



3. LUMINOSITIES, COLORS, AND STAR FORMATION RATES OF 
MERGERS 

3.1. Host Galaxy Luminosities 

To de-convolve the observed pair or merging galaxy lumi- 
nosity function (MGLF) and infer a merger rate as a func- 
tion of galaxy stellar mass, we must first describe the star- 
burst lightcurves. In the simulations, star formation is tracked 
self-consistently, with stellar ages and stellar metallicities de- 
termined directly from the star-forming gas, which itself is 
enriched by star formation. We generally assume that pre- 
existing stars (i.e. initial bulge or disk stars when the gas frac- 
tion is less than unity) were formed with a constant star for- 
mation rate prior to the starting point of the simulation, with 
metal enrichment calculated self-consistently from the ini- 
tial primordial gas. Considering instead exponentially declin- 
ing star formation rates ("r-models") or late-type star forma- 
tion histories fi tted from observations (e.g., Kauffmann et aQ 
'2003b; Ha mmer et al. ''2005) makes relatively little difference, 
and in any case stars formed over the course of the simulations 
usually dominate the observed luminosity. From the predicted 
ages and metallicities of the stars in our simulations, we use 
the stellar population synthesis models of Bruzual & CharloU 
(2003) assuming a Salpetet ( 1955) initial mass function, to 
measure the luminosity in a given band. 

We determine attenuation in the manner described in 
'Hopki ns et all (l2005al l2006bh . Briefly, we calculate the 
column density along ~ 1000 lines-of-sight to each simu- 
lation stellar particle at each time, using the SPH formal- 
ism and multiphase ISM equation of state of the simula- 
tions (Springel & Hernauist 2003a) to determine the density, 
metallicity, and ionization state of gas in the "warm/hot" 
ISM through which the majority of sightlines will pass (i.e. 
neglecting mass collapsed in cold clouds). We adopt a 
gas-to -dust ratio scaled by metallicity (e.g., iBouchet et alJ 
Il985h and normalized to that of the Milky Way, Ab/Nhi = 
8.47 (Z/0.02) X 10-22cm2, and a Small Magellanic Cloud-like 
reddening curve (as suggested by observations for the host 
galaxies of quasars, e.g. Hopkins et al. 2004, Ellison et al. 
2005, although the 280 nm cross-section is decreased by only 
6% if we instead a ssume a M ilky Way-like reddening curve) 
with the form from lPeil (Il992h . to attenuate the intrinsic spec- 
trum of each stellar particle. We do not perform a full ra- 
diative transfer calculation, and therefore do not model scat- 
tering or re-processing of radiation by dust in the infrared. 
However, we compare the column de nsities to quasa r s calcu- 
lated by variations of this method in 'Hopki ns et al.l (|2005b') 
(see their Figures 1, 5, & 6) and find that typical uncertainties 
are a factor ~ 2 - 3 in A^h, generally smaller than the variation 
across different simulations and viewing angles. Furthermore, 
we find below that our results for the visible (attenuated) lu- 
minosities agree well with those expected from the calcula- 
tions in Jonsson et al. ( 2006) which employ a complete Monte 
Carlo radiative transfer model. 

This modeling effectively captures obscuration by gas and 
dust in the "typical" ISM, but does not, however, directly re- 
solve individual molecular clouds or e.g. dense nuclear con- 
centrations of molecular gas, which have important effects on 
the mid and far-IR spectra of the nuclear starbursts of LIRGs 
and ULIRGs. These effects are unimportant in our analysis, 
for several reasons. First, the duration of the starburst is short 
(~ 10** yr) compared to the total time during which a merger 
will be identified 10^ yr), and the fraction of the total stel- 
lar mass formed therein is small. Most of the time the merger 
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Fig. 1. — Ratio of total stellar mass to total A'-band luminosity as 
a function of time in three representative simulations with gas fractions 
/gas = 0.2, 0.4, 1 .0 as labeled (thick lines). For /gas = 1.0, several simulations 
with Vvir = 80-320 kms"' and ^eos = 0.25 — 1 are shown (thin lines). Green 
points show observed ULIRGs (expected to correspond to high-/gas merg- 
ers), with blue points showing the median and interquartile range, from the 
observed merging ULIRG pair sample of Dasvra et al. ( 2006) (r ~ 1.2 Gyr; 
horizontal error roughly illustrates the times corresponding to observed sepa- 
rations), the late-stage ULIRG sample of Tacconi et al. ( 2002 ) (f ~ 1 .5 Gyr), 
and the merger-remnant LIRG/ULIRG sample of Rothberg & Joseph i2004) 
(; ~ 1.9 Gyr). The K-hand mass-to-light ratio depends weakly on gas fraction 
in a systematic manner, but not on other varied merger properties, and agrees 
well with typical observed values. 



will be a weakly-interacting pair or post-merger remnant, so 
corrections to luminosity or mass functions will in any case 
be small, < 10%. Moreover, surveys identifying mergers by 
pair counts or tidal features with which we compare generally 
do not include the heavily obscured starburst stage. 

Second, we explicitly confine ourselves to consider the 
near-IR through near-UV, where the effects of these obscu- 
ration mechanisms are minimized. For typical ULIRGs, cen- 
tral luminosities can be heavily extinguished with A\ ^ 50 
even in /iT-band; but as a consequence the integrated lumi- 
nosity is dominated by stars in more extended regions, which 
are obscured by the relatively quiescent ISM with a much 
lower effective total or "screen" extinction of 0.7 magni- 
tudes (e.g., Lonsdale et al. 2006), similar to that in our sim- 
ulations. Third, re-radiation is relatively unimportant in the 
bands we consider, and at the point of the nuclear starburst, 
even the typical ISM densities are sufficient to render the most 
dense regions optically thick. At such a point, the obscu- 
ration mechanism and the "degree" of optical thickness are 
unimportant. Finally, we have compared with full radiative 
transfer models in the far-IR through UV, including dust heat- 
ing and considering different sub-resolution models for the 
detailed clumping of molecular gas (Narayanan et al. 2006; 
Chakrabarti et al. 2006), and find little change in the total 
luminosity in the bands we consider here, even during peak 
starburst stages. 

3.1.1. K-band Luminosities 

We first consider the /T-band galaxy luminosity in mergers. 
In each simulation, we calculate and plot the total stellar mass 
divided by total /T-band luminosity as a function of time. Fig- 
ure ^ shows the result for several representative cases, some 
with identical virial velocities, but a varying gas fraction, and 
others with the same /gas, but either different Vvi, or qeos- 
There is no clear systematic difference in M/Ly^ at fixed gas 



fraction (we have also varied e.g. initial black hole mass, pres- 
ence or absence of initial bulges, and orbital parameters and 
find the same). However, there is a systematic dependence on 
the gas fraction, so that M/Lk is somewhat lower for high-/gas 
systems which are dominated by younger stellar populations. 

We can estimate this effect quantitatively as a function of 
/gas with a model in which a galaxy of total final stellar mass 
tot forms a fraction /gas of its mass in the merger (i.e. as- 
suming the limit of efficient conversion of gas to stars) with 
mass-to-light ratio (M/LK)new, and the remaining mass com- 
ing from an older population with (M/LK)oid- This gives 
M _ (M/LkXs^ 

Lk~ f +(^-f ^ (^/^k)"" ' 

./gas + U /gas; (M/£K)old 

Fitting to the mean M/Lk in our mergers as a function of /gas 
gives (M/LK)new ~ 0.4 and (M/LK)oid ~ 1 -4, corresponding to 
stellar populations with mean ages ^ 0.5 and 5 Gyr, respec- 
tively, reasonable values for observed mergers. This equation 
gives the mean M/Lk accurate to ^ 10-20%; i.e. comparable 
to the scatter in M/Lk in a given merger and across mergers 
with varied parameters (but similar gas fractions). Note that 
the dependence is fairly weak, as a factor of 2 change in gas 
fraction (e.g. /gas = 0.1 ^ 0.2 or /gas = 0.2 0.4) results in 
only a ^ 20% change in M/Lk, generally smaller than the 
relevant observational uncertainties. 

The ^T-band total mass-to-light ratio is also relatively con- 
stant with time in a merger In part, this owes to the weak 
dependence of M/Lk on stellar population age and the con- 
tinuous star formation throughout the merger For exam- 
ple, given a constant star formation rate (SFR) beginning at 
t = 0, and ignoring older stars (i.e. roughly the /gas = 1 case), 
we expect M/Lk to increase by only ^ 40% from f = 1 Gyr 
to f = 2 Gyr, as in Figure |T] Moreover, as seen in simula- 
tions (Hopkins et al. 2006b; Jonsson et al. 2006) a nd observa- 
tions of local ULIRGs such as Arp 200 dPovon et alJll994t 
Lonsdale et al. 2006), increased densities in the merger give 
rise to both elevated star formation but also enhanced obscu- 
ration, resulting in a flat observed host galaxy luminosity dur- 
ing these times. 

Unfortunately, without a large sample of mergers or a reli- 
able means to determine the exact merger stage or pre-merger 
/gas, comparing our predicted M/Lk as a function of time with 
observations is not possible, but we can check for consistency. 
In Figure ^ we plot M/Lk of observed merging ULIRG 
pairs with masses determined from dynamical measurements 
in Dasvra et al. (2006) at f ^ 1.2 Gyr, corresponding roughly 
to the stage at which the merging nuclear pair separations 
are comparable. Green crosses show individual objects, blue 
points the median and interquartile range, with horizontal er- 
rors heuristically corresponding to times in which a given 
merger would appear in the observed sample. At f ~ 1.5 Gyr 
we show the same, for the late-stage ULIRG sample of 
iTacconi et alJ J2002h . and at t 1.9 Gyr, for the post-merger 
LIRG/ULIRG sub-sample of lRothberg& Joseph (2004). At 
each stage, the typical values and scatter in M/Lk correspond 
well with that expected, especially given the additional scatter 
in observed samples from measurement errors and a range in 
observed /gas (note that the observed samples are all selected 
to be ULIRGs, so we expect they will correspond to the most 
gas-rich mergers we simulate). The later-stage (non-LIRG) 
merger remnant sample of Rothberg & Joseph ( 2004) forms a 
continuum in M/Lk from these values to those typical of old 
ellipticals M/Lk ~ 1-4, presumably as these populations age 
and redden. 



Hopkins et al. 




-22 
















1.0 1.2 1.4 1. 
Time [Gyr] 



-20 -21 -22 -23 -24 -25 -26 

M280 



Fig. 2. — Left: Median observed (attenuated) magnitude at 280nm as 
a function of time during the peak merger stages from three representative 
simulations with Vvir = 80, 160, and 320 kms"' (red, blue, and black, respec- 
tively). Right: Distribution of observed magnitudes (each weighted by the 
amount of time across all sightlines that the galaxy is observed at a given 
M2S0) for each simulation (histograms of corresponding colors), with best-fit 
Gaussians (smooth curves). The distribution of opticalAJV luminosities is 
relatively constant throughout the merger, with enhanced star formation and 
obscuration offsetting one another, and is well-fitted by a Gaussian PDF. 
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Fig. 3. — Left: Mean observed 280 nm magnitude over the course of a 
merger for simulations with Vvir = 80— 320kms"' and gas fractions as labeled 
(circles), as a function of the mass of stars formed during the merger, M, new- 
Vertical errors show the range allowing for a factor of 3 uncertainty in the 
calculated column density. Solid line shows the best-fit power law, dashed 
lines the Icr range of this fit. Right: Dispersion in 280 nm magnitude over the 
course of a merger, for the same simulations (symbols as in left panel). Solid 
line shows the best-fit constant criso and dashed lines the Icr range. The PDF 
for observing a merger at a given luminosity has a constant shape and width, 
with the mean optical/UV luminosity scaling with the mass of stars forming 
in the merger 



3. 1.2. 280 nm Luminosities 

Next, we examine a synthetic UV band centered at 280 nm, 
with a 40 nm w idth, for comparison with the results of 
IWolfetal J (120051) from the GEMS and GOODS surveys. This 
band is normahzed to the solar luminosity in the 280/40 
passband, which is Mq ^w = 6.66 in Vega units or L0.28O = 
2.56 X 10'"W/Hz Here, we use the Stai-burst99 population 
synthesis models JLeitherer et aljfl999) to calculate the UV 
spectra of the stellar particles. Figure 121 shows the outcome 
of this calculation for three representative simulations, with 
/gas = 1 -0, Zgai = 0, ^Eos = 1 -0, and varying Kir, as labeled. 
The median observed M2W as a function of time during the 
merger is shown for each simulation (left panel). Again, the 
effects of star formation and dust obscuration tend to offset 
one another, and the host galaxy luminosity varies within a 
factor < 2, i.e. within 1 magnitude, throughout the merger, 
even though the intrinsic (un-attenuated) UV luminosity can 
rise by a much larger factor. 

We calculate (Figure |2l right panel) the probability of ob- 
serving a given M280, i.e. the total time a given M280 is visi- 
ble integrated over all sightlines and over the duration of the 
merger The results are shown as histograms of the corre- 
sponding color for each merger in the left panel (with the best- 
fit Gaussian). In each case, the PDFs have a well-defined peak 
and narrow (^ 0.6 mag) width in M280, emphasizing that the 
light curves are quite flat. The Gaussian fits allow us to sta- 
tistically characterize the lightcurves of our simulations with 
a mean M280 and rms dispersion about this mean, CT280- 

In Figure |3] we quantify the dependence of the mean ob- 
served 280 nm luminosity on host galaxy properties. We in- 
clude a subset of our simulations with various values of /gas, 
as labeled, and for each, Vyir = 113, 160, and 320kms"' and 
lEOS = 1 -0. For /gas = 1 .0, we also calculate the results for 
Kii- = 80, 226 kms-i, ^eos = 1-0 and for Ki,- = leOkms"', 
Qeos = 0.25, to demonstrate that these variations do not al- 
ter our conclusions. We plot the mean M280 and rms disper- 
sion (T280 for each simulation. The vertical error bars show the 



range resulting from a systematic increase or decrease in our 
estimated column densities (roughly parameterizing the max- 
imal uncertainty for our calculation of A^h from Hopkins et al. 
2005b). We expect the UV luminosity will be dominated by 
young stars, and therefore show the mean M280 as a function 
of the mass of stars formed during the merger, M», new Unsur- 
prisingly, the two are well correlated, and the scaling does not 
depend on gas fraction, virial velocity, or ISM gas equation of 
state. Fitting a relation of the form 

M280 = A/280,() + a log 10 ( , „io . . (2) 

gives a best fit M280.0 = -20.7 ± 0.1, a = -2.01 ±0.13, shown 
as a solid black line. In terms of L28O, 

/ s ().80±0.()5 

L280 =L280,o(y^T3^j , (3) 

where ^280.0 = 5.06 ± 0.45 x IQ^'^Lq 28o- The dispersion (T28O 
has no significant dependence on stellar or black hole mass, 
gas fraction, or other host properties. The solid line (dashed 
lines) shows the best fit (±lcr) mean (T28O = 0.62 ± 0.06 
(0.25 dex in L28o)- There are limited observations to compare 
against, but these values are consistent with the colors and 
inferred star formation properties in e.g. Bell et al. (2006a). 

The scaling of M280 with host properties can be understood 
in terms of a model assuming efficient conversion of gas to 
stars. Given that general trends in star formation rate as a 
function of time are qualitatively similar regardless of the to- 
tal number of stars formed, the intrinsic (un-attenuated) UV 
luminosity (henceforth Lint) should scale proportionally with 
the mass of new stars formed. Lint = (A^/i'28o)~'A^*.new, where 
(M/L280) is an effective mass-to-light ratio depending on the 
star formation history. The timescale of ~ 1 Gyr during which 
the galaxies will be visible as a merger, or alternativel y the fit 
to M/ Lk from our simulations or observations ( Taccon Tet alJ 
l2002h . imply (M/L280) ~ O.O8Mq/L0, 28O, corresponding to 
a population with an age of 0.5 Gyr. 
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However, the observed luminosity does not scale as steeply 
as Lint oc M^.ne.w, owiug t o high er attenuation as the gas mass 
increases. i Jonsson et alJ consider lightcurves of merg- 

ing galaxies using a smaller set of simulations (without black 
holes), but with a full Monte Carlo radiative transfer code, and 
demonstrate that it is reasonable to assume that the newly- 
formed stars and star-forming gas (which produces most of 
the obscuration) are well-mixed, giving an observed luminos- 
ity L28O = iintT"' (1 -e""^) (e.g., Calzetti et al. 1994), where r 
is an effective mean optical depth. From their fits to the at- 
tenuation in bolometric, SDSS M-band, and GALEX NUV lu- 
minosities we estimate r w ro(Mg/lO"M0)" '^, where tq w 
2-3 and Mg w M*_new is the gas mass, which they show can 
be roughly understood in terms of how both the SFR and ob- 
scuration scale with local gas density. Given r > 1, the atten- 
uation goes as L280 oc r"' , and we can combine the scalings to 
obtain the expected relation 

L280 « 5.8-8.7 x10»'Lq, 280 (y^T^j^) , (4) 

in good agreement with the scaling of L28O measured in the 
simulations. That these scalings agree emphasizes that contri- 
butions to the observed L28O from older stellar populations are 
relatively small (< 30%). Furthermore, this suggests that our 
results using just the calculated column densities (i.e. ignor- 
ing scattering processes) would not be significantly changed 
by a more sophisticated treatment of radiative transfer 

Finally, the conversion of gas to stars in mergers is effi- 
cient, i.e. M,^new ~ /gasM*,tot = /gasM* (truc to < 10% ovcr 
/gas ^ 0.2-1 .0), where for simplicity we subsequently denote 
the final, total stellar spheroid mass as = tot- Our cal- 
culations then describe the statistics of A'-band and 280 nm 
luminosities in mergers as a function of /gas and M,, with 
no systematic dependence (quantified in this manner) on the 
other parameters we have varied. 

3.1.3. Other Optical and Near-IR Luminosities 

For future reference, we perform the above set of calcula- 
tions in several different wavebands (280 nm, U, B, V, R, I, J, 
H, K, and SDSS u, g, r, i, z), calculating for each the equiva- 
lent of Equation (|2j, i.e. the best-fit relation between the me- 
dian magnitude in the given band (Mband) and total stellar 
mass formed in the merger (M*_„ew), as well as the disper- 
sion in a given merger about the expected median magnitude 
((Tband)- The results are shown in Table [2 We also consider 
fits to Mband as a function of final black hole mass Mbh or 
total final stellar mass tot instead of M^ ^ew F'or each, we 
show the resulting x^/i^^ which measures the goodness-of-fit. 
Note that a small x^/i^ implies small scatter about the fitted 
relation, but does not imply that there is no systematic depen- 
dence on certain external variables. 

For example, in the optical/UV, the correlation between 
Mband and new or tot is comparable, but as demon- 
strated above, there is a systematic dependence on /gas if 
Mband is quantified as a function of M* tot (brighter opti- 
cal/UV at higher gas fractions). Note that ctband also quan- 
tifies the degree to which it is appropriate to assume that the 
luminosity in the given band over the course of the merger is 
constant (a poor approximation to the bolometric luminosity, 
but a surprisingly good one in the bands we measure). These 
should be of particular use in comparing cosmological simu- 
lations and semi-analytic models with observed merger pop- 
ulations, and enable future observations of merger luminos- 
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Fig. 4. — Color-magnitude relation at ~ 100 random times and sightlines 
during each of several merger simulations with Vvii = 116, 160, 320km s"' 
(diamonds, squares, and open circles, respectively), /gas = 0.2, 0.4, 0.8, 1.0 
(black, blue, green, and red, respectively), Zg-A = 0, ^eos = 1-0. The left panel 
neglects dust attenuation, the right panel includes this effect. In each, solid 
hues show the impact of increasing metallicity (from 0. 1 Zq - Zq ; blue), 
age (from 0.5 — 1 Gyr; green), and column density (from Nn = 0—10^' cm"^; 
red). Filled b lack circles show local ULIRGs from Surace et al. ( 2000) and 
lArribas et iO 120041 . K-corrected following Bell et al i 2005.) . For compar- 
ison, dotted contours (roughly) show the local red seque nce (red) a nd blue 
cloud (blue), with the "valley" separating the two from IBell et ali i |.2005. ) 
(dashed). The color-magnitude distribution of the simulations is similar to 
that observed, with roughly equal contributions from stellar population ef- 
fects and variations in column densities. 



ity functions to easily estimate their contributions to remnant 
spheroid and quasar populations. 

3.2. Colors and Color-Magnitude Relations 

In Figure|3 we plot the (U-V) vs. My color-magnitude re- 
lation for each of several representative simulations with dif- 
ferent virial velocities and gas fractions /gas = 0.2- 1.0, fix- 
ing Zgai and qEOS- For each simulation, we consider ~ 100 
randomly selected points in time at which the merger is "ob- 
served" during the ^ 2 Gyr during which the system shows 
evidence of merging, and at each point in time consider a ran- 
dom sightline (uniformly sampling the unit sphere). In the 
left panel, we neglect the effects of dust attenuation. We illus- 
trate the age/metallicity /obscuration degeneracy with the solid 
lines in the upper right, which demonstrate how a point in the 
color-magnitude space is moved (from the point of origin of 
the three lines) with increasing mean stellar age, metallicity, 
or column density. 

The wide scatter in {U - V) colors in both panels demon- 
strates that the effects of different mean stellar ages and 
metallicities contribute strongly to the scatter in the color- 
magnitude relation. The evolution is difficult to disentangle, 
as increasing age and metallicity in later stages generally red- 
den the merger, but the final merger is associated with a star- 
burst which briefly reduces the mean (luminosity-weighted) 
age. Differential extinction actually marginally decreases the 
color spread, as intrinsically blue periods (i.e. times of peak 
star formation) tend to be heavily dust-reddened, whereas 
intrinsically red periods (e.g. post-merger stages) are only 
weakly affected. Although the scatter in (U - V) may derive 
largely from stellar population effects, this does not imply 
that extinction is unimportant. The mean effective total or 
"screen" reddening is ~ 0.2-0.5 magnitudes and during the 
peak starburst phase this rises to > 1 - 2 mag, similar to that 
observed (e.g., Genzel et al. 1998). (Note, however, that the 
spatially resolved extinction can be much greater in the very 
central regions during these phases.) Extinction also reddens 
objects out of the bright blue (My <-2l,(U-V)<0) region 
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of the color-magnitude diagram (quantitatively, the fraction of 
points in this region drops from ~ 10% to ~ 1% when attenu- 
ation is included). The faint, red region is not especially heav- 
ily populated, which is of interest for optically or UV-selected 
samples, especially at high redshift, which may not be sensi- 
tive to faint red mergers (see e.g. Bell et al. 2005; Wolf et al.l 
I2005L for further discussion). We emphasize though that there 
is no real color-magnitude "relation," and thus any physi- 
cal inferences from the observed color-magnitude distribution 
should be considered with caution. 

We plot for comparison (filled circles) the locations of local 
ULIRGs from Surace et al. (2000) and Arribas et al. (200 J) 
following Bell et aL (.200 5). Quantitative comparison is dif- 
ficult, but the observed and simulated loci are similar, espe- 
cially compared to the blue cloud and red sequence on which 
the progenitor and remnant galaxies lie, respectively (see 
ISpringel et al. 2005a). We have also compared the (280 -V), 
{U - B), (B - R), and (R - K) colors with observed color- 
magnitude distabutions_^^ sample of 
iBell et al.' (200?) and lWolf et all (|2005), and the morpholog- 
ically irregular starburst sample of|de Mello et al. (2006), as 
well as th e GEMS host galaxy s ample of active galactic nuclei 
(AGN) of lSanchez etall (120041) (see e.g. Figure 7 of Bell et al. 
2005), and find similar broad agreement. The substantial scat- 
ter in metallicities during mergers is also comparable to that 
in Rothbera & Joseph (2006), and rough trends in obscuration 
and metallicity with stellar mass or luminosi ty in the star- 
forming stages are simila r to those observe d ("Burste in et al.l 
[1988; Worthev et al. 1992l; lFaber et all 1995t b0rgense nl 19971 
lKuntschneii200Qi) . 

3.3. The Star Formation Rate Distribution in Individual 
Mergers 

Because the star formation rate is far from constant dur- 
ing a merger, sharply peaking during the starburst phase and 
falling off exponentially afterwards, with smaller peaks in ear- 
lier stages, a given merger will not necessarily be observed 
near a single characteristic SFR (as with e.g. the luminosi- 
ties calculated above). We can, however, consider the time 
a given simulation spends in a given interval in SFR (i.e. the 
SFR-dependent "lifetime" of the merger) and use this to quan- 
tify the distribution of SFRs in mergers. In detail, the time- 
behavior of the SFR, M,, is somewhat chaotic, with a system- 
atic dependence on virial velocity, ISM gas equation of state, 
galaxy gas fractions, and merging galaxy mass ratio. How- 
ever, when rescaled in terms of the total mass of new stars 
formed during the merger, M, new, the critical statistical prop- 
erty of the SFR distribution, namely the time spent by a given 
simulation in a given interval in M*, is essentially unaffected 
by these systematics. 

In Figure|5] we determine the integrated time spent by each 
of 100 simulations above a given SFR (from M, ~ 10""^- 
10^ M0 yr"'), as a function of that SFR, M,. We then rescale 
the SFR by dividing out a characteristic SFR, new/?*, i-e. 
the total new stellar mass formed during the merger divided 
by some characteristic starburst timescale. For now, is ar- 
bitrary (we choose f» = 0.65 Gyr, for the reasons below). We 
plot, for each simulation, the time spent above a given dimen- 
sionless SFR (M,/(M*_new/f*))- We show this for the simu- 
lations in bins of log(M, new) from new ~ 10^- lO'-'. In 
each bin, different colored lines correspond to different ini- 
tial conditions in the simulations, various combinations of 
Vvh- = 80-500kms-i, ^eos = 0, 1, Zgai = 0, 2, 3, 6, and 
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Fig. 5. — Integrated time above a given rescaled star formation rate 
M/(M,.„cw/'*), where M»,ncw is the mass of new stars formed during the 
merger and ?, = 0.65 Gyr is fitted. Results are shown for ~ 100 simulations, 
in bins of M*. new as labeled. Different colors show simulations with different 
values of Vyi,-, /gas, 9eos, Zgal- The thick black line in each interval shows the 
fitted analytical approximation of Equation js)- Quantified in this manner, 
the PDF for observing a given merger at some instantaneous SFR is well- 
determined and independent of systematic effects from varied quantities in 
the simulations. 



/gas = 0.2, 0.4, 0.8, 1.0. Although there is significant scatter 
in the time spent above some SFR between different simula- 
tions, we find that once rescaled in this manner, there is no 
systematic dependence on any of the varied quantities. 
We fit the "lifetime" at a given SFR in Figure 13 to an ana- 



lytical function following lHopkins et al 
= f, In(lO) exp(^ 



df 



dlogM, 



(5) 



where new is the total stellar mass formed in the merger, 
is a characteristic timescale, and the normalization In(lO) is 
determined by the integral constraint 



df = / M, 



df 



dlogM, 



^dlog(M,,)=M,^ 



(6) 



This gives a one-parameter function to describe the statistical 
properties of the SFR distribution in mergers. 

Fitting to the simulations, we find a best-fit f* w 0.65 ± 
0.10 Gyr, with no evidence for a strong dependence of f* 
on M* new or other varied parameters. This is compara- 
ble to the star formation times cale im plied by observations 
of starburst galaxies (Kennicut ^ll998h and me an luminosity- 
weighted ages of merger stellar populations dTacconi et alJ 
2002). We show the prediction of this fit as the thick black 
lines in each panel of Figurels] The form of Equation (|5} sug- 
gests an exponentially increasing/decreasing SFR into/out of 
the starburst peak. However, too much detail should not be 
read into the precise SFR as a function of time from these fits, 
as they both average over the entire merger and do not capture 
a systematic dependence such as the tendency of higher-mass 
systems to form a larger fraction of their stars prior to the final 
merger (Robertson et al. 2006a). As an aside, a similar func- 
tional form could be used to describe quiescent SFRs in expo- 
nentially declining "r-models," but in this case would be char- 
acterized by a much lower typical SFR (~ /gas^^disk/lOGyr) 
and higher normalization t In 1 (where t ^ 4 - 8 Gyr is the 
SFR timescale in stable disks; e.g. Li et al. 2005, 2006). 
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4. METHODOLOGY: COMPARING MERGER AND QUASAR 
DISTRIBUTIONS 

4.1. Converting Between Merger, Quasar, and SFR 
Distributions 

If the merging galaxy luminosity function (MGLF) in a 
band i' is given by (j>gaiiL^) = d<i>/dlog(Ly), then 

f dt 

^gdL^) = / . «(M,)dlogM,, (7) 

J dlog(L^,) 

where «(M*) = dndf"' dlog(M,)"' is the merger rate as a func- 
tion of final stellar mass, and df/dlog(Ly) is the time the 
merger is visib le in a logarithmic interval in L,^. Because, 
as shown in § 13.11 the merging galaxy luminosity is quite 
flat in time, and the timescale over which merging galaxies 
are identifiable as such 2 Gyr) does not depend strongly 
on the galaxy masses over the range of interest, we can use 
df/dlog(L„) = P{L^ M ^,)f m erge, whcrc P{Li, |M,) is the prob- 
ability (calculated in § 13. 1> of being observed in a logarith- 
mic interval in over the merger timescale (fmeige)- The de- 
convolution of the MGLF, (piL^), with the PDF for observing 
a merger at a given luminosity P{Li, |M,) directly yields the 
merger mass function for a given luminosity function. 

Having considered the observed properties of mergers with 
some final spheroid mass M,, we can use the scalings be- 
tween black hole mass and host galaxy velocity dispersion 
or stellar spheroid mass to relate each merger to an associated 
black hole. The relation between final (post-merger) black 
hole mass and total stellar mass in our simulations is 

Mbh^O.OOIMh,, (8) 

jDi Matteo et al .''2Q05': 'R obertson et al.l '2006b'). with an ap- 
proximately lognormal scatter of ~ 0.3 dex, essentiall y iden- 
tical t o observational estimates from e.g. (Marconi &Hun3 
J2003h : iHaring & Rixl (12004) (accounting for the difference 
between virial and stellar mass). Knowing Mbh(A^*), we can 
convert n(M,) to niM-^n), the birthrate of quasars of a given 
relic mass in galaxy mergers, 

n(MBH) = J P{Msii\M,)ri(M,)d\ogM,. (9) 

We can then express the observed QLF in terms of «'(Mbh) 
and the quasar lifetime as a function of luminosity (see Equa- 
tion IIIOI ). identical to the expression in Equation Q but with 
M, replaced by Mbh- 

Given a model for the quasar lifetime df/dlog(Li,), 
then, the QLF is determined. We generally consider two 
cases: first, a "feast or famine" or "light-bulb" model 
for the quasar lifetime, following what has generally 
been adopted in previous works (e.g., 'Small & Blandford' 
1^92; Kauffmann & Haehnelt 2000; Haiman & Menou 
2000; Wvithe&Loeb 2003; IVolonteri et alJ I2003t 
Haiman. Ouataert. & Bower ,20Q4.) . In one typical sce- 
nario of this type, the quasar turns on in a merger, accretes 
at constant Eddington ratio m, and then either "turns off" 
or ex ponentially decays in luminosity (Haiman &Loe3 
Il998h . The lightcurve is an exponential, L cx exp(±f/fQ), 
and the lifetime as a function of luminosity is trivial, 
df/dlog(L) = In 10 fg = constant (for L < Lpeak), where tg is 
a free parameter. Assuming an even simpler description in 
which quasars "turn on" at fixed luminosity L = Lpeak for a 
fixed tQ, produces nearly identical results. Generally, tg is 
either adopted from observations (with the loose constraint 
tQ 10^- 10^ yr; see Martini 2004 for a review) or assumed 



to be the Salpeter (1964) time ts = 4.2 x 10^ yr (e-folding time 
at m = 1 with radiative efficiency er = 0.1), but we allow it to 
vary to provide the best fit to the observed QLF. 

In Hopkins et al. (2005a,b, 2006a,b), we compare these 
idealized scenarios to the light curves and lifetimes derived 
from our simulations and find that they provide a poor repre- 
sentation at any luminosity. The simulation lightcurves are 
complex, with multiple phases of activity, variability, and 
feedback-regulated decay, but their statistical nature can b e 
described by simple forms ( Hopk ins et al.l l2005b 2006tf). 
The key feature of the quasar lifetime not captured by more 
idealized models is that it increases with decreasing luminos- 
ity; i.e. a given quasar spends more time at luminosities be- 
low its peak than at its peak luminosity. In Hopkins et ak] 
c) we use our simulations to calculate the differen- 
tial quasar lifetime, i.e. the time spent in a given logarithmic 
luminosity interval, and find it is well-fitted by a Schechter 
function, 

^-.(^)""exp[-L/L*]. (10) 

Here, Lg is proportional to the peak quasar luminosity Lpeak ~ 
LEddiMjiH), tQ is a fixed constant {tQ = l^^fs for L*q = rjL^e^^), 
and a ^ 0.5 is weakly dependent on peak luminosity and is 
determined by the nearly scale-invariant "blowout" of gas as 
exponentially growing feedback from black hole growth heats 
the gas rapidly and it can no longer cool efficiently in a dy- 
namical time ( Hopkins et al. 2006a). The lifetime in Equa- 
tion ( I10> is entirely determined by our simulations, and when 
quantified as a function of Lpeak in this manner, the quasar life- 
time shows no systematic dependence on host galaxy proper- 
ties, merger parameters, initial black hole masses, ISM and 
gas equations of sta te and star formation m odels, or other var- 
ied parameters ('see lHopkins et al.l ( I2006ti) for a detailed dis- 
cussion and calculation). 

Having quantified the luminosity dependent quasar or 
merger "lifetime" in a manner independent of the detailed 
initial conditions above, we can map between merger and 
quasar populations without requiring any cosmological priors 
on their distributions. Thus while we do not calculate the dis- 
tribution of host galaxy properties a priori but rather model in- 
dividual mergers in order to achieve the resolution necessary 
to model processes related to star formation and black hole 
growth, by quantifying our scalings in terms of e.g. Mbh and 
tot, our methodology is independent of the vast majority of 
host galaxy properties. These differences manifest in e.g. dif- 
ferent final stellar or black hole masses, but do not change any 
scalings expressed in terms of stellar or black hole mass. It is 
straightforward to reverse this procedure, for example start- 
ing with the QLF and quasar lifetime to determine liiM^n), 
converting to ri{M^), and using this to predict the MGLF. We 
can also compare with other quantities, e.g. the SERF, by con- 
volving n(M*) with the SFR-dependent Ufetime f(M* |M*). 

4.2. Systematic Uncertainties: Calibrating Merger 
Timescales/Selection Effects and Gas Fractions 

Before explicitly comparing merger, quasar, and SFR dis- 
tributions as a function of redshift, we consider systematic 
uncertainties introduced in the various observed samples by 
e.g. selection effects and cosmic variance. For quasar popula- 
tions, these should not be significant sources of uncertainty, as 
we generally employ hard X-ray samples which are complete 
to luminosities well below those of interest in our analysis. 
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especially at the moderate red shifts which we consider in de- 
tail (see e.g. lUeda et an i2003). and have been corrected for 
obscuring columns with Nh S 10^^ cm"^. We have also de- 
signed our formulation for comparing e.g. quasar and merger 
populations to be robust with respect to the wide set of param- 
eters varied in our simulations, as discussed in §|3] 

However, for a certain "true" merger rate as a function of 
mass, selection effects will influence the observed merger 
mass and luminosity functions. Of course, there will be some 
mass/luminosity completeness limit, but it is straightforward 
to account for this. More problematic is the fact that differ- 
ent definitions of mergers (i.e. choice of maximum separation 
in pair-selected samples or strength of tidal features in mor- 
phologically identified samples) will, for the same distribu- 
tion of merger rates, select different sets of "ongoing" merg- 
ers. Fortunately, we find in §|3that merger light curves are 
quite flat as a function of time, with no strong dependence 
of observed total luminosity (in the near-UV through near- 
IR bands we consider) on merger stage. This means that the 
observed distribution of merger luminosities and masses will 
not be biased by the various sample definitions, only the total 
number of objects identified as mergers (i.e. vertical normal- 
ization of merger luminosity functions). This uncertainty is 
entirely contained within the parameter fmeige defined above; 
i.e. the total time over which a system would be identified as 
a merger in an observed sample. In principle, our simulations 
can be used to determine this quantity from first principles 
for a given observational sample, convolving predicted galaxy 
images with the appropriate observational response (includ- 
ing e.g. surface brightness dimming, sky and detector noise) 
and selection criteria, to calculate the probability along each 
sightline that the object will be identified observationally. 

Such an investigation is beyond the scope of the current pa- 
per, but is clearly important for detailed estimates of merger 
rates as a function of mass and redshift for a uniformly se- 
lected merging galaxy sample at different redshifts. Some 
preliminary results can be obtained by tracking the path of our 
simulated mergers in Gini-M20 space, a non-parametric mor- 
phological classification scheme in w hich spirals, ellip ticals, 
and mergers occupy distinct regions fL otz et aU mrn. This 
quantitative analysis implies that the galaxy interaction may 
be classified as a merger somewhere between ^ 0.5-3.0Gyr, 
depending on the viewing angle, the orbit, the progenitor disk 
orientations, and the dark matter halo concentration. The me- 
dian merger time is w 2 .0 Gyr and it is worth noting that a sim- 
ilar estimate obt ains from dyn amical time considerations and 
observations (P atton et al. !'2000). In what follows, we simply 
adopt this as the (constant) fmeige for our comparison with ob- 
served merger luminosity functions. This is no doubt a rough 
approximation, but we can test whether variations in fmeige will 
have a significant impact on our conclusions. 

Our mapping between e.g. MGLF and QLF is also not en- 
tirely independent of our simulation parameters, in the sense 
that a given QLF defines a unique merger mass distribution, 
but does not constrain the initial (pre-merger) gas fractions 
of the progenitor galaxies. For a given mass tot, the char- 
acteristic luminosity at which the merger is observed can be 
affected by /gas, at least for optical/UV samples. Furthermore, 
we expect that the distribution of gas fractions will be drawn 
from some PDF, which may depend on galaxy mass. We 
do not model this (instead adopting a mean ,/gas ^ given 
redshift), primarily because it introduces additional, poorly 
constrained parameters and because it is unclear to what ex- 
tent this correlation is driven by more massive systems having 



undergone previous major mergers (i.e. their pre-merger gas 
fractions may have been uniform). However, the effect is less 
important than might be expected, for two reasons. 

First, the scaling with mass is not steep (i.e. a factor ^ 
10^- lO-' in M, yields a factor ^2-3 in /c,as, similar to the 
existing u ncertainty ). Second, our modeling and the obser- 
vations of lXuetalJ (2004); Wolfetal. (2005); Bund v et alJ 
J2005ah imply that the merger mass and luminosity functions 
do not have steep low mass/faint-end slopes, meaning that it 
is a good approximation to adopt the typical mean (/gas) of 
galaxies near the MGLF "break". In detail, we have consid- 
ered our calculations a dopting the f^as as a func tion of late- 
type galaxy mass from lRoberts & Havnesl(ll994h . and find it 
makes little difference to our predictions (generally modify- 
ing them within existing ^ Icr uncertainties). Of course, to the 
extent that /gas evolves with redshift (see below), this consid- 
eration is important in disentangling how much of that evolu- 
tion owes to a change in the characteristic masses of merging 
galaxies or to gas depletion through quiescent star formation. 

Even if the vertical and horizontal normalizations of the 
MGLF (equivalently, fm ei-ge and /gas) were completely unde- 
termined, we show in S IsVll that the MGLF and QLF shapes 
are not trivially related, and this alone provides a strong test 
of models of quasar lightcurves and triggering. Moreover, 
these systematics (contained in /gas and fmerse) are reasonably 
well-constrained as described above, and the systematic de- 
pendencies in our comparisons are not especially strong over 
the range in which they are constrained. 

We directly infer the sensitivity of our predictions and com- 
parisons between e.g. the QLF and observed MGLF at a given 
redshift to these systematics. In §|5]we use the formalism de- 
scribed in § 14. H and our simulations to map the observed QLF 
to a predicted MGLF, and find that it agrees we ll with ob- 
served MGLFs at z ~ 0.2 in the g-b and (Xu et al. H2004ft and 
z - 0.7 at 280 nm f Wolf et al.ll2005& . for expected values of 
fmerge ~ 0.5 - 3 Gyr and /gas ~ 0. 1 - 0.6. We can quantify this 
formally: for any given "mean" f^erge and /gas at some red- 
shift, the QLF and our methodology define a predicted MGLF, 
and we can determine the probability (in a sense) that the 
observed MGLF measures this predicted MGLF. 

In Figure|6l we plot this for the observed MGLFs at z ~ 0.2 
and z ^ 0.7, where probability contours (essentially measur- 
ing the agreement between observed and predicted MGLFs) 
are shown in fmeige and /gas. In the /T-band (z ~ 0.2), this com- 
parison is essentially independent of /gas (as expected from 
the weak dependence of /T-band luminosity on /gas in Equa- 
tionQ, and depends only weakly on fmerge in the range fmerge ^ 
l-2Gyr At 280nm, the dependence on /gas is stronger, as 
expected since L28O scales with M,.new ^ /gasA/*, but our re- 
sults are still essentially unchanged for reasonable values of 
typical z ^ 1 quasar-producing mergers with gas fractions 
/gas ~ 0.2-0.6 and the broad range of fmeige ^ 0.5 -3 Gyr. 
Even considering the effects on both MGLFs simultaneously 
in the right panel of Figure |6l (which plots the joint prob- 
ability contours from the /T-band and 280 nm MGLFs), our 
comparison is unaffected over the range fmerge ^ 1 -2Gyr and 
/gas ~ 0.2-0.6 (these numbers should not be taken too liter- 
ally, since the samples derive from different redshifts). 

In future samples which determine the merging galaxy lu- 
minosity function by uniform, automated selection criteria 
(see, e.g., Lotz et al. 2006), our comparisons can be made 
more precise by inferring fmerge a priori from the simulations 
and directly measuring /gas in optical samples. The key point 
from the above, however, is that our theoretical factor ~ 2 
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Fi g. 6. — Pro bability contours in liost galaxy gas fraction /i.as and time the merger is visible in the observed sample ^mcree for the MGLF in AT-band at z < 0.2 
fleft: IXu et aTi2 004) and at 280 nm at z = 0.7 (middle; Wolf et al. 2005). Contours are shown at 10, 25, 50, 75, 90% as labeled, and quantify the degree to which 
systematics will change our results. Right panel shows the previous two contour sets (dotted) and the joint (combined) probability contours of both luminosity 
functions (solid). Allowing these to vary freely yields physically reasonable values, and comparison of the two demonstrates that the dominant systematic 
uncertainties are different in near-IR and optical/UV MGLFs. Systematics from selection effects and theoretical uncertainties should not change our results over 
the reasonable ranges of fmergc and /gas, and will generally be smaller than the effects of cosmic variance. 



uncertainty in fmerge (which stems from not modeling the full 
sample selection function and e.g. cosmological orbit distri- 
butions) is not important for our conclusions, as any value 
of fmerge in this range still produces a good fit to observed 
QLFs (i.e. gives a self-consistent mapping between QLF and 
MGLF). Likewise, variations in merging galaxy gas fraction 
do not change our results. 

Moreover, cosmic variance is likely to overwhelm the sys- 
tematic normalization issues above for the small volumes 
probed in observational measurements of merging galaxy 
statistics. We discuss this further in § 16.21 below, but note 
that observations with the deep imaging required to identify 
merging syste ms have generally been limited to small fields. 
ISomerville et al. ( 2004b) and l^Tf e t al. ( 2005) consider this, 
and estimate that cosmic variance introduces a systematic fac- 
tor of ~ 1 .5 - 2 in the normalization of the merger luminosity 
functions. Therefore, it is not warranted to employ a more 
detailed treatment of merger sample selection effects (fmerse) 
or /gas in our modeling, which considers the relation between 
e.g. the MGLF and QLF, if these are measured in different 
fields and either (more seriously the MGLF) is significantly 
affected by cosmic variance, and in any case the introduced 
systematic uncertainties are comparable to current measure- 
ment errors in observed MGLFs. 

5. COMPARING THE MERGING GALAXY, QUASAR, AND STAR 
FORMATION RATE DISTRIBUTIONS 

5.1. The Observed MGLF and QLF 

5.1.1. The K-Band Merger Luminosity Function 

We examine first the observed /iT-band pair luminosity func- 
tion of Xu et al. (2004), at low redshift, z < 0.1. The pair lu- 
minosity function is determined fro m the matc h ed 2M ASS- 
2dFGRS 45,289-galaxy sample of ICole et al.' ('200!|, and 
agrees well with p reviou s estimates of the local pair fr action 
(e.g., Zeof&Koo 1989; Burklev et al. 1994; Carlbers et alj 
[1994; .Yee & Eningson 1995; Pattonetal. 1997. 2000) and 
the fi-band luminosity function of paired galaxies (e.g.. 
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2003), but accounts for effective pair volume corrections in- 
stead of treating each member singly with its own Vmax (since 
both objects must be identified in a pair to be included in 
such a sample). Pairs are defined to be within a projected 
separation < 20 /i"' kpc, with velocity difference (where both 
redshifts are measured) < 500km s"'. Further, the sample 
is restricted to major mergers, with A'-band magnitude dif- 
ferences less than 1 mag (i.e. within a factor ~ 3 in mass), 
corresponding well to our modeling because these mergers 
are most likely to trigger starburst and quasar activity and 
to be visible as peculiar/interacting galaxies for com parison 
with m orphologically sele cted MGLFs ( Mihos & Hernauisll 
fT9943 lWalker etaljll996t IConselice et al 2003) . We con- 
sider both their binned data and best-fit Schechter function, 
with slope a = -0.30 ±0.56, turnover M/^. ,^ = -23.32 ±0.25, 
and normalization log 0, = -3 .92 ± 0. 1 3 (converted to our cos- 
mology). The distinction between the observed mass func- 
tion and a constant merger fraction is only marginally signif- 
icant; moreover, a number of studies measure a cumulative 
merger fraction. Therefore, we also consider our predictions 
for a constant merger fraction as a function of mass (given 
the local stellar mass function from iBell et a l.' f2003")) with 

the observed low -redshift valu e ^ 2.5 ± 1% tXu et al 20041 

l^ndy et al. 2005a b; Bell et a l.l2006bh . 

Figure shows (upper left panel) the predicted QLF 
from the /T-band low redshift MGLF of Xu et al. (2004). 
Solid ci rcles in dicate the observed z < 0.2 hard X-ray 
QLF of lUeda et al. (2003J. We rescale the hai'd X-ray 
QLF to a b olometric Q LF for ease of comparison as in 
lHoDk ins."Richards. & Hernquisll (|2006'), using a model of 
the intrinsic quasar continuum SED following Marconi et 
(^004), based on optical through hard X-ray observ ations 

tg., Elvis et al. 1994; Georee et al. 1998; Vanden Berk et ajj 
01; Perola et al. 2002; Telfer et al. 2002; Ueda et al. 2001 
gnali et al.ii2003i) . with a hard X-ray reflection component 
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Fig. 7.— Left: Predicted QLF from the observed ^T-band z < 0.2 rupperi lXuet alJ2004ft and 280 nm z, a 0.7 aoweri lWoif et alJ20O.'5;) MGLFs, adopting our 
full model for quasar lifetimes (solid black line; blue shading shows Icr range based on errors in the observed MGLF), or a "light-bulb" model for the quasar 
lifetime with tQ varied to produce the best-fit expected QLF (short dashed; green hatch ed region shows Icr range; dotted lines show prediction for tQ a factor of 10 
larger or smaller). Circles show the observed hard X-ray quasar luminosity function of lUeda et ali i2QQ^), rescaled to bolometric luminosity with the corrections 
from|Marconi et al. i2004). The predicted QLF from a constant merger fraction 2.5 ± 1% at z < 0.2 is also shown (dot-dashed line uses our full model, 
long-dashed line a "light-bulb" model). Middle: Predicted AT-band z < 0.2 and 280 nm z~OJ MGLFs from the QLF of Ueda et al. ( 2003), lines and colors show 
the same as the left panels. Observations from Xu et al. (2004) from the 2MASS sample are shown as black squares in the upper panel, and a constant merger 
fraction (2.5 it 1%; given the mass functions of Bell et al. 2003) as circles. In the lower panel, observations from Wolf et alj t200.5i) are shown as black squares 
(GOODS), blue diamonds (GEMS), and red circles (GEMS depth on GOODS area). Right: Inferred fi(M»), the gas-rich merger rate (spheroid birthrate) with a 
given stellar mass, from the corresponding (AT-band, upper; 280 nm, lower) observed MGLFs (Icr range, blue horizontal hatched; instead using a constant merger 
fraction shown as blue vertical hatched), and the corresponding QLF aXz< 0.2, z ~ 0.7, respectively (Icr range, red shaded). Lines with circles and squares show 
the fitted ri(M*) from the QLF used in Hookins et al. 1 2006a b), respectively. In all cases, our modeling demonstrates that the observed QLFs and MGLFs are 
self-consistent and can be used to predict one another, with all mergers producing bright quasars and all bright quasars initially triggered in mergers. Idealized 
models of the quasar lightcurve give a very misleading relation between MGLF and QLF, predicting a QLF discrepant by > 6 cr. 



gen erated by the PEXRAV model JMagdziarz & Zdziarskil 
1199 5^. This includes the observed dependence of the opti- 
cal to X-ray luminosity ratio on luminosity (e.g., Wilkes et alJ 
^94; Green et al. 1995; VignaH et al. 2003; Risaliti & El^ 
llpOS; Strateva et al. 2005), but for our purposes adopting 
the mean observed quasar SED from lRichards et all J2006b l) 
yields a similar result (see e.g. Hopkins et al. 2005d, 2006b 
for a discussion of the impact of different bolometric correc- 
tions on our predictions). For simplicity, we consider only the 
hard X-ray (2-10 keV) Q LF, where atten uation is usually neg- 
ligible (and note the .Ueda et alJ i2003h QLF is corrected for 
columns A^// < 10^^ cm^). 

The black solid line in Figure0shows the prediction using 
the quasar lifetimes fitted from our simulations (described in 
§ 14. U . with the blue shaded region giving the la range al- 
lowed based on the errors in the observed MGLF (dot-dashed 
line is the result adopting instead a constant merger fraction). 
We show the results for a characteristic fmerge = 2 Gyr, corre- 
sponding roughly to the time when the galaxies are within 
^ 50kpc of one another in our simulations before merging 
(appropriate for the maximal projected separation 20/1"' « 
30kpc of the observed sample), and /gas = 0.15, corresponding 
to the typical gas fractions of galaxies around the mass of the 
observed MGLF peak ( Roberts & Havnes 1994). However, as 
shown in § 14.21 our results are relatively insensitive to these 
choices. Briefly, if we allow /gas and fmerge to vary, a formal 
best-fit is obtained for /gas = 0.25^0 23 and fmerge = l-7^a9 Gyr 
- in other words, the quality of the fit is nearly unchanged for 



the entire reasonable range of /gas and fmerge, and these spe- 
cific values should not be taken too literally. The short dashed 
line shows the prediction adoptin g the idealized "light-bulb" 
quasar model described above (§ 14. 1> . with again the green 
horizontal cross-hatched range showing the la range from 
the observed MGLF (long-dashed line with vertical hatching 
shows the result for a constant merger fraction). The predic- 
tion shown in this case is the best-fit allowing fmerge, /gas, and 
the quasar lifetime Iq to vary freely, giving fg w 6.3 x 10* yr 
for fmerge = 2 Gyr (only the ratio fg /fmerge is constrained) and 
/gas = 0.25. The dotted lines show the prediction for tg larger 
or smaller by an order of magnitude. 

In our model, the faint end of the observed QLF is dom- 
inated by quasars with large peak luminosities, but seen in 
fainter evolutionary states during and after the mergers in 
which they form. Therefore, the faint end of the predicted 
QLF is tightly constrained, as it depends on the MGLF near 
its turnover, where it is most well-determined, and there is 
little difference between adopting the observed MGLF or a 
constant merger fraction. With instead a one-to-one cor- 
respondence between peak and observed quasar luminosity 
(dashed lines), the faint-end QLF depends directly on the 
faint-end MGLF, and is thus poorly constrained. The agree- 
ment between the observed QLF and our prediction is good, 
X^/i^ ~ 3.41/8 = 0.43 and it requires no fine-tuning of any 
parameters. However, even allowing both horizontal and 
vertical normalizations to vary freely, the shapes alone of 
the quasar and MGLFs rule out a "light bulb" or exponen- 
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tial lightcurve model at greater than 99.9% confidence, with 
the best-fit y} jv « 29.1/7 = 4.15 (x^/z^ « 6.45 adopting a 
constant merger fraction). Such a fit furthermore requires a 
quasar lifetime tgr^Gx 10^ yr, at the low extr eme of the range 
indicated by observations (e.g.. Martini 2004). 

In any model which reproduces the observed Mbh oc re- 
lation, «(Mbh) must have approximately the same shape as the 
merger rate as a function of M», itself given approximately by 
the merger mass function. This shape is observationally well- 
defined, with a rapid decline above and flattening or decrease 
below the break luminosity. In Hopkins et al. (2005c), we ar- 
gue that this is the shape of «(Mbh) (equivalently «(Lpeak)) 
implied by the combination of our model for quasar life- 
times with the observed QLF. However, in idealized models 
of the quasar lightcurve, «(Mbh) must have the same shape as 
the observed QLF, which cannot be self-consistently resolved 
with the observed MGLFs. It is true that the faint-end slope of 
the MGLF is poorly constrained, with a = 0.30 ± 0.56; how- 
ever, to be compatible with the faint end slope of the QLF in a 
light-bulb model (i.e. to steepen the dashed line to match the 
solid line in Figure requires a 4-a change in a (6(7 relative 
to a constant merger fraction), which seems unlikely, given 
the large quoted error which includes systematic effects. 

Next, we invert this comparison and employ the observed 
QLF to estimate the MGLF. Figure \7\ (upper middle panel) 
shows this, using the z < 0.2 hard X-ray Ueda et al. (20031) 
QLF shown in the upper left panel to determine n(MBH), 
and correspondingly n(M,) and (pgaiiLx), the ^T-band z < 0.2 
MGLF. Again, the blue range shows the prediction and la 
range (from errors in the observed QLF) using our full model 
of quasar lifetimes, with /gas = 0.15, fmerge = 2Gyr Dashed 
line and green hatched range shows the same for the "light- 
bulb" quasar lifetime model, allowing /gas, fmeige, Iq to be 
fit as in the upper left panel (producing identical fits), and 
the dotted lines show the prediction for an order of magni- 
tude larger or smaller tg. These can be compared to the black 
squares, which show the binned MGLF from Xu et al. ( 2004J) 
(open circles show the constant merger fraction estimate). As 
expected from the upper left panel, the agreement using our 
full model of quasar lifetimes is good jv = 0.35), and the 
agreement with an idealized model of the quasar lifetime is 
poor (xV^ = 3.82). 

The constraints from the observed QLF on the MGLF are 
significantly weaker than the constraints on the QLF from the 
observed MGLF. This is not because the QLF is poorly con- 
strained relative to the MGLF (in fact, the opposite is true). 
Rather, it is because a given interval in observed quasar lumi- 
nosity has significant contributions from quasars with a wide 
range of peak luminosities Lpeak ^ L (and a correspondingly 
wide range in host galaxy stellar mass), in various stages of 
evolution. Thus, there are significant degeneracies in predict- 
ing the relative contributions from different black hole or host 
galaxy stellar masses to the observed QLF based only on its 
observed shape. This is increasingly true at fainter merger lu- 
minosities, and therefore we can place only weak constraints 
on the faint-end ri(MsH) distribution (as discussed in e.g. Lidz 
et al. 2006). The weak constraints at the high-L^: end, how- 
ever, owe to larger statistical errors in the QLF. The wide 
range in the allowed faint-L^ end of the MGLF (demonstrated 
by the similar predictions given the observed MGLF or a con- 
stant merger fraction) is reassuring, in the sense that at low 
luminosities, selection effects become important in efforts to 
identify a sample of merging galaxies, but these will not sig- 
nificantly change our result. 



As an intermediate stage in these calculations, we have de- 
rived the implied «(M,); i.e. the merger rate or birthrate of 
spheroids with stellar mass, directly related to the black hole 
"triggering" rate. We show these n(M*) distributions in Fig- 
ure (upper right panel), where the blue horizontal shaded 
range sho ws the la range implied by the observed MGLF of 
IXu et alJ (1^04) (vertically shaded range implied by a con- 
stant merger fraction), and the red range shows the la range 
impHed by the observed QLF of Ueda et al. (2003) at z < 0.2. 
Again, it is clear that the QLF provides poor constraints 
on the low-mass behavior of «'(M*). For comparison, we 
plot the fitted «(M,) distributions of Hopkins et al. (2006a b), 
which were fitted to the combination of optical, soft X-ray, 
and hard X-ray Q LFs (Mivaiietal. 2001; Ueda et al. 20031 
'Croom et al ."2004; Richards et al. 2005; Hasinger et al. 20051 
La Franca et al. 2005). These may rise to lower masses before 
turning over than is implied by the MGLF, but as we have em- 
phasized above (and shown in e.g. Hopkins et al. 2006b) this 
produces an essentially identical prediction for the properties 
of the quasar and remnant red galaxy populations. 

5.1.2. The 280nm Merger Luminosity Function 

We repeat t he analysis in § 15.1.11 for the observed 280 nm 
UV MGLF of lWolf etal.l ll2005) from the GEMS and GOODS 
surveys at z « 0.7. In detail, the galaxies are selected from 
the overlapping area of COMBO-17 and GEMS, with photo- 
metric redshifts 0.65 < z < 0.75 and observed R <24 (corre- 
sponding approximately to a rest-frame U -band). The sample 
consists of visually identified morphologically selected pecu- 
liar/interacting galaxies, and as a consequence there is consid- 
erable observational uncertainty in the faint end of the MGLF, 
where tidal features are especially difficult to detect given sur- 
face brightness dimming at higher redshifts. However, as we 
demonstrate above in the case of the A'-band MGLF, our mod- 
eling of the QLF is not sensitive to this faint-end behavior 

Figure \7\ shows this comparison in the lower panels, in 
the same manner as the upper panels. The lower left panel 
shows the predicted QLF using our full model for the quasar 
lightcurve (solid line), and instead for an idealized "light- 
bulb" model for the quasar lightcurve (dashed line), with un- 
certainties as in the upper panels, and dotted lines show the 
prediction for a fixed quasar lifetime an order of magnitude 
larger o r smaller. We com pare to the observed hard X-ray 
QLF of lUedaetal.1 (l2003l) . this time at z = 0.7 (black cir- 
cles). Whether we use the GEMS or GOODS data makes 
little difference in the prediction of our model, so we show 
the result using the GOODS luminosity function, which has 
a larger number of faint merging galaxies, maximizing the 
ability of the "light bulb" model to match the data. Again, 
it is clear that our model provides a good fit, while the 
"light bulb" model is ruled out at greater than 90% confi- 
dence jv = 0.159, 2.41 respectively). The shape of the 
QLF alone provides this qualitative constraint on the quasar 
lightcurve model, but the vertical and horizontal normaliza- 
tions also agree well within the reasonable expected ranges 

fmeige ~ 1 - 3 Gyr and /gas ^ 0.2-0.6. As above, fmerge, /gas, 

and fg are allowed to vary freely in the "light bulb" model, 
and although the fit is poor, the best-fit values are /gas = 0.3 
and fg = 5.8 x 10^ yr for fmerge = 2 Gyr. A ~ 4ct change in the 
faint-end slope of the MGLF would be required for the "light 
bulb" model to agree with the QLF, even for the ma ximal es- 
timate of the number of faint mergers in Wolf et al. (2005). 

The lower middle panel of Figure 171 shows the 280 nm 
MGLF inferred from the z = 0.7 QLF of lUeda et all (l2003l) . 
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in the same style as the upper mid dle panel. The in ferred 
MGLFs can be compared to those of IWolf et al.l (120051) from 
the GOODS survey, GEMS survey, and a combination of 
GEMS depth on the GOODS area. The differences in the ob- 
served MGLFs at the faint end demonstrate the considerable 
observational uncertainty here, but the QLF is not sensitive 
to this faint-end behavior This is not the case for the ideal- 
ized "light bulb" model, as demonstrated by the narrow range 
of the green hatched region even at low luminosities. The 
constraints are weaker, but the agreement between our model 
predictions and the observed MGLF is good (y^ jv = 0.31, 
compared to 'X^ lv = 2.15 for the light-bulb model). 

Finally, we show the n(M^) distribution implied by the 
MGLF of Wolf et al. (2005) (lower right panel) as the blue 
shaded range. The range implied by the QLF at z = 0.7 is in- 
di cated as the re d shaded range, with the fits to «(M*) used 
in lHopkins et al.l (l200 6a b). The implied n(M») distributions 
at z = 0.7 and z < 0.2 from 280 nm and /T-band observations, 
respectively, can be compared to infer the evolution of merger 
rates and quasar birthrates as a function of mass, but we cau- 
tion against taking this comparison too far The samples use 
various selection criteria, are from different surveys, are not 
in the same wavebands, and have different dependence on gas 
fraction. Therefore, the potential systematic effects may be 
quite different between the two. 




Fig. 8. — Merger mass function at z ^ 0.7, pred i cted fr om our model- 
ing witli the observed 280 nm MGLF of Wolf et alj 12005) (colored), and 
directly calculated from near-IR and optical observations in Bundy et aD 
i2005a) (black). Points converted from the 280 nm observations in Figure 171 
are shown as cyan stars (GOODS), blue diamonds (GEMS), and red circles 
(GEMS depth on GOODS area), and the ~ Icr allowed mass function from 
our modeling as the yellow shaded region. Black squares show the binned 
mass function from Bundv et al. 1 2005a), black dotted line shows the best-fit 
Schechter function. The agreement implies that our simulations reliably map 
between MGLFs and merger stellar or black hole mass functions. 



5.1.3. Merger Mass Functions 

We repeat our a nalysis using t he morphologically separated 
mass functions of iBundv et al.l ( 12005 a^b.) from GOODS and 
DEEP2. Stellar masses are determined by fitting to optical 
spectroscopy and near-IR photometry, and we adopt the mass 
functions for objects visually classified as mergers (following 
iBrinchmann et alJll998l) with ACS morphologies. This gives 
a binned merger stellar mass function determined in each of 
the intervals z = 0.2-0.55, z = 0.55-0.8, and z = 0.8-1.4, 
spanning 10^- lO'^M© (see Figure |8] below for a de- 

tailed comparison with the z = 0.55-0.8 mass function). In 
each case, the observed merger mass functio n implies a QLF 
in good agreement with the observations of Ueda et alJ ( l2003h 
(xVi^ wO.85, 0.52, 0.98 for z = 0.2-0.55, 0.55-0.8, 0.8- 
1.4, respectively), and the implied mass function from the 
QLF, albeit less well-constrained, is consistent with the ob- 
served mass functions. These observations also consistently 
rule out a "light bulb" or similarly idealized model for the 
quasar lightcurve at > lOa. 

5.2. Comparison of Merger Mass and Luminosity Functions 

While the agreement between quasar and merger statis- 
tics is encouraging, we now test our modeling of the 
MGLF directly by comparing the observed 280 nm MGLF of 
Wolf et al. (2005) to the measured merger mass funct ions of 
Bundv et al. (2Q05,a b). We consider the Bundv et aL (|2005i) 
mass functions determined from the GOODS field, with red- 
shifts 0.55 < z < 0.8, similar to the 0.65 < z < 0.75 range con- 
sidered in Wolf et al. ( 2005). In both cases, optical or near-IR 
morphologies are used to visually select merging systems, so 
the ambiguity in selection effects is minimized. 

Figure|8]shows the inferred merger mass function using our 
simulated P(L28o|-W*) to convert 4>{L2U)) directly to (j>(M^). 
The observations of Wolf et al. (2005), converted to a mass 
function, are shown by the symbols indicated in the cap- 
tion. Alternativel y, adopting the ^ la range in allowed fit- 
ted MGLFs from IWolf et all (120051) and converting this to a 



merger mass function yields the shaded yellow range. The 
range shown can also be thought of as a range in gas fraction, 
with /gas ~ 0.15-0.30. The observed merger mass function 
of Bundy et al. (2005a) is shown as the black squares, with 
the black dotted line showing the best-fit Schechter function. 

The merger mass function implied by the MGLF and our 
modeling agrees well with that observed, suggesting that we 
are capturing the critical dependence of observed merger lu- 
minosities on host properties. Coupled with observations of a 
similar relation between black hole and spheroid mass at these 
redshifts (e.g., Peng et al. 2006), this observationally supports 
multiple implicit levels of our modeling. 

5.3. The Distribution of Star Formation Rates 

Having determined «(M*) from the QLF or MGLF, we can 
convolve instead with the time each merger spends in some 
interval in SFR, (Equation (5)), to determine the observed 
SERF in mergers. Figure |9] shows the results given n(M*) 
(specifically «(M* new)) determined from th e 280 nm observed 
peculiar/interacting luminosity function of Wolf et al.l (120051) 
at z = 0.7 from GEMS and GOODS (soHd line with Icr blue 
shaded range). We compare this to the obs erved SFRF of th e 
same merging galaxy sample, determined in lBell et al.l(l2005l) . 
Note that while Bell et al. ( 2005) determine the SFRF and rel- 
ative number density as a function of morphology, they do 
not determine the absolute SFRF normalization (i.e. proper 
ef fective volume correcti on), and so we adopt that inferred 
bv lLe Floc'h et alJ ( 12005ft at z = 0.6-0.8 from Spitzer in the 
Chandra Deep Field South. Red points are observations with 
no correction for incompleteness below ^ lOMoyr"', 
while black points have been corrected for incompleteness at 
low SFR based on the c omparison of GEM S and GOODS lu- 
minosity functions from' Wolf et alJ ( 12005 ). 

The agreement is good, implying that we are properly mod- 
eling the relation between the optical/UV luminosity of merg- 
ers and at least the statistics of the induced star formation as 
a function of time. This is not trivial, as there is no one-to- 
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Fig. 9. — Distribution of star formation rates in m ergers at z = 0.7, i nferred 
from our modeling and tlie z = 0.7, 280 nm MGLF ofWolfet al.1 12005) (solid 
line). Shaded range shows Itr range based on observational errors in the 
MGLF. Circles show the observed SFRF of Bell et al. 1 2005) at the same red- 
shift (black, corrected for completeness; red, uncorrected), normalized with 
the observed IR luminosity functions of Le Floc'h et al. 1 2005). Red line 
shows the predicted SFRF owing to major mergers from the semi-analytic 
models of Somerville et al. i2001), blue Hnes include minor mergers. SoHd 
hnes adopt ?mcrsc ~ lOr^yn ~Gyr, dotted lines rjy„ ~ 10* yr as the maxi- 
mum time since the last merger The agreement imphes that our modeling 
reliably predicts the merger-driven SFR distribution from a given MGLF or 
QLF. Note the increasing importance of minor mergers at low SFR. 



one correspondence between the 280 nm luminosity and SFR 
(as optical/UV lightcurves are flat while SFRs change by or- 
ders of magnitude). Also note that, because both the 280 nm 
magnitude and SFR distributions are functions of new (as 
opposed to M*_tot), the relation between the two is indepen- 
dent of M*,tot and /gas. 

It is also of interest to compare our prediction for the distri- 
bution of mergers and SFRFs to that from semi-analytic mod- 
els. In general, this is difficult, as the events of interest are not 
mergers of dark matter halos, but involve the luminous galax- 
ies themselves, in which the galaxies are of comparable mass 
and have a large supply of cold (rotationally supported) gas. 
Furthermore, the exact requirements for triggering starbursts 
and AGN depend on, for example, the pressurization of ISM 
gas, star formation recipes, and the distribution of orbital pa- 
rameters in the mergers (e.g., Hernquist & Mihos 1995). We 
have specifically chosen to quantify the properties of mergers 
in terms of e.g. Lpeak, Af*,tot, and „ew, in a manner which 
suppresses these dependencies, in order to relate merger, star- 
burst, and quasar populations independent of these detailed 
(and difficult to model) cosmological distributions. It is there- 
fore outside the scope of this paper to consider a detailed com- 
parison of merger rates as a function of these conditions as 
implied by e.g. our semi-empirical modeling, semi-analytic 
models, and cosmological simulations. However, we briefly 
consider a comparison between our predicted merger-driven 
SFRF and that predicted by semi-analytic models in Figure|9l 

We show as the red solid line the SFRF predicted owing to 
major mergers (mass ratios less than ~3:1), and as the blue 
solid line the SFRF owing to "minor" mergers (mass ratios 
less than ~10: 1), calculated at z = 0.7 from the semi-analytic 
models of Somerville et al. (2001, 2004a). We consider a 
galaxy to be "merging" if it has undergone a merger within a 
time ~ lOfdyn '^Gyr (blue dotted line shows the minor- merger 
prediction for a merger time ^ t^yn ^ 10** yr). There is broad 
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Fig. 10. — Contours of the Icr allowed Schechter function parameters M«, 
(assuming a = -0.5) for the MGLF in A'-band at z < 0.2 (left) and at 
280nmatz = 0.7 (right), assuming different values for the gas fraction /gas 
and selection efficiency (time during which the merger will be identified as 
such) fracigc, as labeled. Points with error bars show the measurements at the 
corresponding frequency and redshift of Xu et al. ( 2004) from the 2MASS- 
2dFGRS survey in ^-band and Wolf et al. 1 2005) from GOODS (square) and 
GEMS (circle) at 280 nm (observations rescaled assuming a = —0.5). The 
QLF can significantly constrain the allowed MGLF parameters (p* and M, , 
modulo systematic uncertainty in the merger timescale fmerge (i.e. sample se- 
lection effects) and typical merging galaxy gas fraction /gas. 



agreement, suggesting that the rate of gas-rich mergers in the 
standard CDM cosmology is appropriate to that implied by 
our modeling and the observations. However, it is only with 
the inclusion of some "minor" mergers that the semi-analytic 
model predictions match the observations and our predictions 
from the QLF and MGLF, suggesting that our implied merger 
rates as a function of final, total stellar mass may have contri- 
butions from a wide range of mass ratios ^10:1 to 1:1 (which 
we implicitly account for by expressing the merger rate in 
terms of the final mass and peak quasar luminosity), increas- 
ingly dominated by major mergers at the high-mass end. This 
comparison, though, is sensitive to the estimate of the ob- 
servable timescal es in the semi-analytic model (but see e.g., 
iBell et alJl2006bL for comparison with these merger rate esti- 
mates). 

6. THE EVOLUTION OF MERGER LUMINOSITY FUNCTIONS AND 
STAR FORMATION RATE DISTRIBUTIONS 

6.1. Merger Luminosity Functions 

Now, we use the observed QLF to predict the MGLF at 
redshifts where measurements do not yet exist. We have seen 
that the observed QLF is insensitive to the number of low- 
mass mergers. If we assume that the MGLF is well-fitted by 
a Schechter function, 

^(M) = 0, io"-4(«-*^-)("+i) exp { - 10°4(^-^*) } , 



(11) 



this means that the QLF has almost no power to constrain 
a. However, if we assume a specific a, we can then fit the 
observed QLF to determine the allowed range of (f)^ and M» . 

FigurefTolshows this determination of (p^ andM*, assuming 
a constant a = -0.5, for both the A'-band z < 0.2 an d 280 nm 
z = 0.7 MGLFs, compared to the observations of iXu et alJ 
ll200l and lWolf etaiT ( l200 5l). The observed MGLF 0* and 
M* have been re-fit assuming the same a = -0.5. For each 
wavelength, we show the contours enclosing the Icr range of 
(/) ^,, M ^, for a given gas fraction and fmerge- As expected from 
§ 14.21 the theoretical uncertainty in fmerge^ /gas is generally 
unimportant relative to the observational factor ~ 2 uncertain- 
ties. In principle, the number of observed mergers scales with 
the selection efficiency fmerge and their luminosities with /gas, 
such that they independently control the prediction of and 
M*, respectively. If there were no errors in observations of the 
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relevant luminosity functions, comparison of M» alone could 
fix the distribution of /gas. However, even small errors in the 
observed luminosity functions make 0, and strongly de- 
generate parameters, and the choice of a introduces another, 
albeit well-known, degeneracy. 

Given these degeneracies, we focus on the more robust pre- 
diction of the luminosity density. The full luminosity density, 
Pl = 0* L, r(a + 2), technically includes the order unity cor- 
rection from the integrated faint-end contribution, r(a + 2), 
which our modeling has little power to constrain from the ob- 
served QLF. However, the quantity 0* is well-determined, 
and in Figure^Jwe show our predictions for this combination 
(i.e. neglecting the r(a+2) correction, or equivalently assum- 
ing a = -1 or a = 0, although the corrections from a are small 
for the reasonable range -1 .5 <a< 0.5). Most observational 
estimates of pi have been in the B-band, so we show the pre- 
dicted luminosity density in /T-band (left) and B-band (right), 
as a function of redshift fr om z = 0-6, based on the combined 
iMivaii et al. (2001); Ueda et alJ dfO OS): Croomet al. (2004); 
l^chards et al. (2005); Hasinger et al. (2005); La Franca et alj 
12005) QLFs. The B-band (440 nm) calculation of host galaxy 
lumi nosities is i dentical to our calculation at 280 nm (see 
also ' Hopkins et aLi20061\ for a full calculation); the lumi- 
nosities and scaling with M,_new are nearly identical (see 
Table Q modulo a nearly mass-independent Mb - M280 ~ 
0.9 (reflecting the mean young stellar population age ^ 
0. 5 Gyr). W e comp are with the z < 0.2 MGLF measurement 
of IXu et akl (120041). the merger IR and B-band luminosity 
densities of Conselice et al. (2003, 200|) at mean redshifts 
z = 0.88, 1.68, 2.22, 2.71, the peculiar/interacting galaxy 
mass functions of Bundv et al. (2005a b) at z = 0.2-0.55, 
z = 0.55-0.8, and z = 0.8-1.4, th e z < .1 optical pair lu- 
minosity function of 'Toledo et al.' (flOOp), and B-band p e- 
culiar luminosity densities from Brinchman n et al.l (Il998h at 
z = 0.2 - 0.5, z = 0.5 - 0.75, and z = 0.75 - 1 .0. 

In order to predict pi at each redshift from the QLF, we 
must adopt a typical gas fraction /gas(z). In Figure fTTI the 
dotted line is for a constant /gas = 0.3, where this value is 
chosen because it provides a best fit to the cumulative ob- 
servations plotted (considering the fit to the full MGLF where 
possible). The solid line assumes an exponentially declining 
gas fraction with cosmic time, /gas = 10"'/'" where fn is the 
Hubble time, such that /gas = 1 at early times and /gas = 0.1, 
similar to the Milky Way, at present. This gives an e-folding 
time for /gas of / ln( 1 0) « 6 Gyr, similar to that expected for 
quiescent spirals following a Schmidt-type star formation law 
jkennicutt 1998; Rownd & Young 1999; Martin & Kennicujj 
1200 L .Springel & Hernquist ,2003a: Li et al.. ■2005.. .200^1 
Comparing with the evolution of the distribution of cold (rota- 
tionally supported) gas mass fractions in galaxy mergers (with 
total baryonicmass above lO'^M©) in the semi-analytic mod- 
els of Som erviile et a l. (2004a), we find that these cosmolog- 
ical models predict a similar evolution in the mean gas frac- 
tion. The yellow shaded range shows the Icr allowed range 
inferred from the QLF, for the evolving /gas case (uncertain- 
ties are similar for constant /gas). 

Although it is clear that the observations do not strongly 
distinguish between the two cases, and we caution that we 
have assumed a constant selection efficiency fmerge despite 
different observed sample selection criteria and different red- 
shifts, there is a significantly better fit to the combination of 
luminosity density and MGLF observations in the case of the 
exponentially declining gas fraction ix^ /v = 0.4 compared to 



/v = 2.1). This distinction is derived from the B-band and 
280 nm data; the /T-band is relatively insensitive to /gas and 
does not significantly constrain this evolution. That the ex- 
ponential decline in gas fraction is preferred by the data is 
expected, and the observations are not accurate enough to sig- 
nificantly constrain the timescale or detailed functional form 
of the decrease in gas fraction with time, but this demonstrates 
the key qualitative behavior, that merging galaxies were more 
gas rich in the past and that these measurements can constrain 
that evolution through the evolution in M»(z). 

It is clear from the figure that our predictions for the lu- 
minosity density are narrowly constrained, generally within 
a factor '--^ 2 at all redshifts, with the errors dominated by 
uncertainties in the observed QLF (and thus potentially im- 
proved by future, high redshift complete samples which can 
constrain the QLF break luminosity at a range of redshifts) 
as opposed to fitting degeneracies. The variation owing to 
different gas fractions is small in the /T-band, as expected, 
but enters approximately linearly into the B-band luminosity 
densities. Improved constraints on the /T-band MGLF mea- 
surements and extension of B-band MGLF measurements to 
higher redshifts can provide a strong test of whether the con- 
sistent relation between merger and quasar populations this 
modeling demonstrates at z ^ 3 remains true at high redshifts. 

6.2. The Merger-Driven Star Formation Rate Distribution 

Following § I5.3l but beginning from the observed QLF in- 
stead of the MGLF, we can predict the SFRF as a function of 
redshift. Figure^jshows the results of this calculation at sev- 
eral redshifts z = 0-3 (as labeled), from the fit to the combined 
Mivaii et al. (2001); Ueda et al. (2003); Croom et al. (200% 
[Richards et al. (2005); Hasinger et al. (2005); La Franc a et all 
12005) QLFs. The dashed line shows the prediction from 
the MGLF (as opposed to the QLF) at z = 0.7 from Fig- 
ure |9] Points a t z = 0.7 sh ow the observations dBell et alJ 
2005; Le Floc'h etani2005l) as in Figure |9l We adopt our 
rough estimate /gas = ~ exp(-f/6Gyr). The evolution 

in the SFRF is predominantly a reflection of the luminosity 
evolution of the QLF, as the break in the QLF evolves to- 
wards larger luminosities with redshift, implying that typi- 
cal corresponding black hole and host galaxy masses must 
evolve similarly. Our assumption that the typical gas frac- 
tions rise with redshift also increases the SFR for a given to- 
tal (final) mass, resulting in a further factor ~ 3 increase in 
new from z = 0.2-2. Comparing e.g. the specific star for- 
mation rate distributions as a function of stellar mass, age, 
and redshift from our simulations and observed in Bauer et alJ 
(12005); Feulneretal ( 2005) ; Perez-Gonzalez et al. ( 2005<) ; 
iPapovich et alJ ( l2006h yields a similar self-consistent "down- 
sizing." 

Although this calculation suggests a form for e.g. the 
ULIRG and LIRG IR luminosity functions, we caution that 
we do not model re-radiation of absorbed light by dust and are 
not yet in a position to predict these properties. However, a 
rough calculation shows that we expect ULlRGs counts to be 
consistent with our modeling (Hopkins et al. 2006b). We de- 
fer a more thorough calculation of the IR spectrum, including 
dust heating, re-radiation and line emission, to future work. 

We can integrate the SFRF calculated at each redshift to de- 
termined the SFR density in mergers as a function of redshift, 
as shown in Figure [O] The solid black line is our prediction 
assuming /gas = 10"'/'" as estimated in § 16.11 with shaded la 
range based on the errors in the observed QLF and degenera- 
cies in fitting. We also show the predictions for a constant 
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Fig. 1 1 . — Predicted luminosity density contributed by merging galaxies in A'-band (left) and B-band (right), as a function of redshift. Solid line assumes a 
gas fraction /gas = 10"'^'", dotted line assumes a constant /gas = 0.3 at all z- Yellow shaded range shows the Icr allo wed range inf erred from the QLF, for the 
evolving fgas case (errors are similar for constant /gas). Observations are shown from Xu et al. i2004) (red triangle), IBundv et alJ l2005a b) (blue diamonds), 
l&inselice et al. (2003, 2005) (circles; HDF-N in cyan, HDF-S in green), Toledo et al. ( 1999) (red squares), and Brinchmann et al. ( 1998) (blue stars). The 
observed mass/luminosity density in mergers and its evolution is consistent with that required if all bright quasars are triggered in mergers and the converse, that 
all mergers trigger bright quasars. The quasar luminosity function and our modeling can be used to predict the merger mass and luminosity density at all redshifts 
where the QLF break is reasonably constrained. 
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Fig. 12. — Infen'ed SFRF in mergers from the QLF as a function of redshift (black lines) at t he redshifts shown. Th e shaded area shows the Icr range from 
errors in the observed QLF and uncertainties in fitting to it. Observations of Bell et alj 12003):*Le Floc'h et al.l •20()3) are shown at z = 0.7 in the manner of 
Figure 1^ and the dashed line in each panel shows the prediction from the 280 nm MGLF of Wolf et al. 1 2005) from Figure|51at z = 0.7 for comparison. The 
dashed line and sohd line atz = 0.7 are slightly different because the former is predicted from the observed merger luminosity function at that redshift, the latter 
from the observed QLF. The QLF alone allows us to reasonably constrain the distribution of star formation rates in mergers, over a wide range of redshifts where 
observations are not available. 
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Fig. 13. — Star formation rate density in mergers, inferred from tlie quasar luminosity function assuming an average gas fraction /gas = W'/'" = exp(-r/6Gyr) 
(solid line). The shaded area shows the Icr range allowed (ow ing to degeneracies in fitting to the observed QLF). Note that adopting the steeper faint-end slope 
estimates from the deeper GOODS subsample of IWolf et alj |2005) systematically increases our prediction by a factor ~ 50% (at the upper range of the la- 
uncertainty). The dot-dashed Une and dashed lines show the result assuming a constant gas fraction /gas = 1.0 and /gas = 0.1, respe ctively (for clarity, errors ar e 
not shown, but they have identical form). Dotted lines show the estimated total (extinction-corrected) star formation rate density from Hopkins & Beacom ( 200f^ 
Observations of the SFR density in mergers at various redshifts are shown from Bell et al. ( 2005 ) (stai), Brinchmann et al. ( 1 998 ) (squares), Perez-Gonzalez et aj 
|5005) (x's; from the estimate of the ULIRG contribution to the SFR density), and Menanteau et al. i200L .2006) (filled circles). The /gas = 1 line can also be 
thought of as the total rate at which stellar mass is "moved" or "generated" onto the red sequence by quasar-producing gas-rich mergers. The QLF alone can be 
used as an independent constraint on the merger-driven SFR density of the Universe; the SFR density in mergers must rise from ~ 1% of the total SFR density 
atz = to ~ 15-25% atz = 1 and ~ 30-50% atz = 2. 
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/gas =1-0 (dot-dashed) or /gas = 0.1 (dashed) (errors similar 
to the case shown). The solid line effectively interpolates be- 
tween a high-gas fraction era at z > 2 to present characteristic 
low gas fractions /gas 0.1 at z = 0. The predicted evolution 
in p(SFR) is understood from the QLF. The time-averaged 
growth rate of black holes is ji/erC^, where ji is the quasar 
luminosity density, so if all bright quasars (i.e. most of jl at 
moderate/high redshift) come from spheroid-forming mergers 

with M,.new ~ /gasA^* ~ /gasA^Be/^ ~ 0.001) WC CXpCCt 

p(SFR)-/gas;L/Me.c2. 

For comparison, we show the total (integrated over all 
morpholog ical types), extinction-corr ected SFR density es- 
timated in Hopkins & Beacom' ('2006') (similar to e.g. Cole 
et al. 2001, Hopkins 2004 or predictions of cosmologi- 
cal simulations, e.g. Springel & Hernquist 2003b, Hern- 
quist & Springel 2003) as the dotted line. Observations 
of the SF R density in mer ging/peculiar system s are plot- 
ted, from Bell et alJ J2005I). iBrinc hmann et al.' ('1998'), and 
iMenanteau et all ll2001[ |2006h. T he estimate of the SFR 
density in ULlRGs from Perez-Go nzalez et a l. (2005) is also 
shown, but we caution that while these objects are associated 
with mergers at low redshift, dusty, more concentrated disks 
at high redshift may contribute to this population. 

The agreement between the observations of the SFR density 
in merging galaxies as a function of redshift and our model- 
ing is good, implying that the QLF can be used as an inde- 
pendent, albeit indirect constraint on the star formation his- 
tory in mergers and the contribution of merger-driven star- 
bursts to present stellar populations. Systematic uncertain- 
ties in the observations are still a concern; e.g. the observa- 
tions of Brinchmann et al. ( 1998) may be incomplete owing 
to dust extinction; see Menanteau et al. (2006). More impor- 
tant, cosmic variance in likely to overwhelm the systematic 
normalization issues above for the small volumes probed in 
the observations shown. The difficulty of identifying mergers 
without deep imaging has naturally limited the effective vol- 
um e of most such m easurements. Somerville et al. (2004b) 
and lWolf et alJ ( l2005 r) discuss this for several of the fields we 
have considered (e.g., GOODS, GEMS, HDF), and for the 
specific application to different morphological classifications, 
and they estimate a factor of ^ 1.5-2 uncertainty in their es- 
timates from GOODS. Lacking a sufficiently large volume to 
overcome cosmic variance, we should ideally compare mea- 
surements from the same fields. That is, regardless of cosmic 
variance in a particular field, the observed QLF, MGLF, and 
merging galaxy SFRFs in that field should be consistent. In 
fact, if we limit our comparisons in this manner, we find very 
good agreement between the observations and our predictions 
(see e.g. §|5|i. The problem here is that again, existing fields 
with deep space-based imaging are too small to contain rare, 
luminous quasars. In general, our predictions from the QLF 
may be more representative of the cosmic mean, as they agree 
with the observed QLFs from large volume surveys such as 
the SDSS (e.g., Richards et alJ2005Ll2()()6l . 

Despite systematic uncertainties, the trend in /3(SFR) ap- 
pears robust, and we can make several comparisons between 
the SFR d ensity in mergers a n d the total SFR density esti- 
mated by Hopkins & BeacomI ( 1200 6). Uncertainties in this 
comparison are dominated by the factor ^2-3 scatter in 
observational estimates of the total SFR density, rather than 
uncertainties in our predictions. Regardless, we confirm the 
results of various observational estimates of merger frac- 
tions and their evolution with redshift ( Le Fevre et a i] l200(l 
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namely that at z < 1, star formation trig- 
gered in mergers contributes only a small fraction to the SFR 
density of the Universe, but this increases near the era of peak 
merger and quasar activity at z > 2. Specifically, from Fig- 
ure|Olwe predict that only ^ 1% of the SFR density at z = 
is driven by mergers, and even assuming a maximal /gas = 1.0 
in all mergers increases this only to ~ 10-20% (setting a 
strong upper limit). The predicted contribution from merg- 
ers rises to ^ 15-25% by z = 1, with an upper limit '^35% 
(/gas = !)■ At z = 2, the contribution from mergers increases 

+10c 



h of the cumulative SFR density, and 35:^"% at z = 3. 



to 40!^' 

Above these redshifts, the best-fit prediction approximately 
preserves these ratios, but the uncertainties become large. 
Thus, spiral galaxie s (and incre asingly irregular galaxies at 
higher redshift, e.g. ICross et alJ (|^04); Daddi et al. (2004); 
ISomerville et al.l (l200 4a^') dominate the SFR density at z < 2, 
even assuming the most optimistic /gas =1.0 pure gaseous 
mergers, with the contribution from mergers rising to a com- 
parable but not dominant ~ 30-40% at higher redshifts. 

These predictions do depend systematically on the faint-end 
slope of the MGLF, or equivalently the low-mass slope of the 
merger mass function, ^fe have adopted typical measured val- 
ues from Xu et al. (2004); Bu ndv et all (2005a b), but if the 
low-mass incompleteness is large and we instead adopt the 
maximally steep slope from Wolf et al. (2005), this increases 
our predictions by ^ 50%, enhancing the merger-driven con- 
tribution to p(SFR) to ^ 30% at z = 1 and 60% at z = 2. 
Moreover, /gas is somewhat uncertain, but it enters approxi- 
mately linearly in Figure^] so our prediction can be rescaled 
with improved observational constraints on /gas. 

By convolving our predicted rate of spheroid formation as 
a function of redshift with the characteristic gas fraction as 
a function of redshift, we can estimate a "mean" /gas for 
spheroid-producing mergers leading to spheroids of a given 
z = mass. The redshift-dependent rat e of spheroid for- 
mation is given in iHopkins et alJ ( l2006bh . but is essentially 
determined by the merger rates as a function of total, final 
mass at each redshift. For example, the /gas = 1 line in Fig- 
ure^lcan also be thought of as the total rate at which stellar 
mass is "moved" from the blue sequence to the red sequence 
by quasar-producing, gas-rich mergers. Because the implied 
merger activity peaks at z > 2 corresponding to the peak of 
quasar activity, assuming /gas = 10"'/'" gives a nearly con- 
stant characteristic mean gas fraction ~ 60% for essentially 
all spheroid masses of interest (M* > 10^ Mq). The old ages 
of spheroids coupled with the expected increase in disk gas 
fractions at high redshift means that it is indeed a good ap- 
proximation to assume that most spheroids are formed in quite 
gas-rich mergers. This is supported by observations, as e.g. 
Hoekstra et al. (2005) find from stellar population analysis 
that elliptical galaxies must have formed much of their stellar 
mass around the time of the elliptical-producing event. This 
does not, of course, prohibit subsequent spheroid-spheroid 
mergers, which will not trigger star formation or quasar ac- 
tivity. 

7. QUASAR HOST GALAXY LUMINOSITY FUNCTIONS 

Given the luminosity function of merging galaxies which 
produce quasars, it is relatively simple to convert this into 
the expected luminosity function of quasar "host galax- 
ies" (HGLF), but we must first define what we mean by 
quasar "hosts." If we include e.g. starbursts with buried. 
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low-luminosity X-ray AGN then we will recover a similar 
merger luminosity function to that we calculate from the X- 
ray QLF. However, most efforts to measure the distribution 
of AGN host galaxy luminosities (e.g ., Bahcall et al. 1993 
McLure et al. 1999; Falomo et al. 2001; Hamilton et al. 20o3 
Jahnke & W isotzki 2003; Dunlo£et al. 2003; Flovd et all 



Vand en Berk et al. 200^h ave considered hosts of 
bright, optical, broad-line quasars. In our modeling, this 
phase of quasar activity is associated with the "blowout," i.e. 
the final stages of black hole growth, when the surrounding 
gas is expelled and heated and the black hole is briefly ren- 
dered a bright optical source before accretion shut s down. 
We consider this in detail in iHopkins et al.l (1200 6b*) (Equa- 
tion [26]) and derive the useful mean relation from our simu- 
lations 



j-qso 

^B, obs 



:7.9 



1 



-/ 



/gas \10^Mq/ Lpeak 



0.2 



(12) 



We can use this to estimate the optical luminosity function 
of quasar hosts by taking our already determined optical/UV 
merging galaxy luminosity functions and adding the selection 
criteria appropriate to these galaxies being identified as the 
hosts of bright optical quasars. 

We show in § 13.11 that the median merger optical/UV lu- 
minosity is roughly constant in time, but the quasar lumi- 
nosity can vary by orders of magnitude and spends much 
of its time at luminosities well below Lpeak ^ LEdd(A^BH)- 
Consequently, the observed quasar luminosity is uncorrected 
with the host galaxy l uminosity. This is discussed further in 
iHoDkins et alM OOSc. •2006b') in the context of the observed 
lack of correlation between quasar luminosity and black 
hole mass (e.g.. Ho 2002; H eckman et alJ 120041: iHao et all 
120051) . but similarly explains the observed lack of a corre- 
lation between nucl ear and host galaxy optical luminosity 
JBahcall et al."199'A 'McLure et al. ffl999t LTahnke & WisotzkJ 
12003; Haoetal. 2005; Vanden B eA et al.1 l2006h . At the 
brightest luminosities above the break in the QLF, this is no 
longer strictly true, as black holes do increasingly tend to be 
at a high Eddington ratio and near their peak luminosity, and 
a correlation between nuclear and host luminosity is expected 
(although at these luminosities measuring host properties is 
most difficult). 

Using our modeling of the quasar and galaxy luminosities 
during a merger, we can predict the joint distribution in opti- 
cal quasar and host galaxy luminosity. In Figure we gen- 
erate a mock sample of ^ 2000 optical "AGN" and plot their 
quasar and host galaxy Z?-band magnitudes. The procedure 
by which the m ock dis tribution is generated is described in 
detail in Hopk ins et al.l (^006b), but briefly, we consider in- 
tervals of log(MBH) from the distribution of relic black hole 
masses n(MBH) at z = 0.5 with 10*'Mbh < Lpeak < lO'^M©. 
For each, we consider the simulations with the nearest Mbh, 
and calculate the PDF for their being observed as optical 
AGN (defined here by Mg. qso < -16 and an obscuring column 
density Nh < 10^^ cm"^) with a given AGN and host galaxy 
B-band luminosity (both calculated including attenuation us- 
ing the methodology in Hopkins et al. (2006b), but adopting 
the empirical fit to column density distributions as a function 
of luminosity from Ueda et al. (2003) yields very similar re- 
sults). From this distribution, we randomly calculate 2000 
points in Mb, qso, Mb, gai, which are shown in the figure. We 
assume /gas = 0.2, but this only changes the normalization 
of Mb, gal in this distribution. For comparison, the dashed 
line in the figure shows Lb, gai = Lb, qso and the solid line 




B, qso 



Fig. 14. — Predicted optical B-band quasar and host galaxy magnitudes 
for a mock sample of 2000 optical "AGN" with Lpeak > 10' L0, at red- 

' O.lLEdd), solid 



shift z = 0.5. Dashed line shows Lg 



■al 



^B, qso 



(igso 



line Lg =0.1 Lb, qso (^QSO ~ ^Edd)- Quasar and host luminosities are un- 
corrected except at the brightest luminosities where systems are generally 
near-Eddington, as observed (e.g.. .Hao et alj2Q05i) . 




Fig. 15. — Predicted B-band HGLF of optical quasars at z = 0.2. Dashed 
line shows the full best-fit B-band MGLF (assuming a = 0.5), solid lines show 
the HGLF requiring an observed B-band quasar luminosity above a given 
fraction of the host galaxy luminosity (left panel, as labeled) or requiring an 
observed B-band quasar magnitude above some fixed limit (right panel, as 
labeled). Points show the observational estimate of iHamilton et al. 1 200^) 
(requiring My, qso < —23, open circle s) at z ^ — 0.4 , rescaled to typical B- 
band host luminosities following Bahcall et al.VlQQV) and Vanden Berk et ^ 
(2006). Although these quasars are not necessarily predicted to be observed 
in an obvious merger stage, their host luminosity function is a subset of the 
MGLF, selected at a different merger stage. 



shows Lb. gal = 0.1 Lb. qso- These approximately correspond 
to Eddington ratios of L w 0.1 Lead and L « LEdd, respec- 
tively. The distribution demonstrates the lack of correlation 
between quasar and host galaxy luminosity, and agrees well 
with various observational estimates (e.g., Bahcall et al. 199 
McLure et al. 1999; Jahnke & Wisotzki 2003; Sanchez et al 
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2004). Specifically, compare Figure 13 of Vanden Berk et al. 
2006, who find a nearly identical distribution considering host 
galaxy-AGN spectral decomposition of SDSS AGN complete 
to Lg. gal ~ 0.1 Lg. qso (where the SDSS g-band is generally 
equivalent to B-band for our purposes). 

For a given Mbh (or correlated merging galaxy luminos- 
ity), we can calculate from Equation il2\ and our scaling of 
the quasar lifetime as a function of instantaneous and peak lu- 
minosities the time spent above some observed quasar B-band 
luminosity limit Lb, qso and/or some limit in the ratio of quasar 
to host galaxy optical luminosity Lb, qso/Ls. gai; i.e. the time 
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during which we estimate the merging galaxy can be identi- 
fied as a traditional optical or broad-line AGN host galaxy. 
Figure[2]shows this calculation assuming different selection 
criteria in determining the HGLF. In the left panel, we show 
the MGLF at z = 0.2, determined from the z = 0.2 hard X-ray 
QLF of Ueda et al. ( 2003) as the dashed line. For clarity, we 
do not show the uncertainty in this prediction, but adopt the 
best-fit assuming a = -0.5 and /gas = 0.2; the uncertainty in 
the plotted MGLF can be seen in the upper and lower middle 
panels of Figure As solid lines, we then show the HGLF, 
i.e. the MGLF given some selection criteria for the quasar and 
host galaxy B-band luminosity, namely Lb, q^o > O.ILb. gai. 
Lb, qso > Lb, gai, and Lb, qso > IOLb, gai, as labeled. In the right 
panel, we show the same, but instead imposing the Z?-band 
quasar magnitude limits Mb, qso < -19, Mg, qso < -21, and 
Mb. qso < -23, as labeled. In each, the observationally esti- 
mated optical quasar HGLF of Hamilton et al. (2002) is plot- 
ted, which is composed of quasars primarily of optical mag- 
nitude Mb < -23 (crosses with error bars; open circles show 
the inferred luminosity function applying a strict My < -23 
magnitude cut). There is some uncertainty in the absolute 
normahzation of the host galaxy luminosities, discussed in 
iHamilton et all (12002^: therefore we rescale their V-band host 
galaxy luminosities to the B-bandusing the ca li brations of the 
overlapping samples in Bahcall et al. (1997); M cLure et al.l 
( fT999) . which gives a peak in the HGLF similar to that ob- 
se rved in e.g . Bahcall et al. ( 1997) and the much larger sample 
of lVanden Berk et al.. (^06) at Mb ~ -22. 

In any case, when we require that a bright optical quasar 
(Mb < -23) be observable as such, we find good agreement 
between our predicted HGLF and the observed distribution of 
host galaxy luminosities. There are several uncertainties in 
this comparison: for example, the normalization and absolute 
host magnitudes of the observed sample may be biased if the 
quasar light is not perfectly subtracted, our prediction only 
loosely estimates the full distribution of dust obscuration as a 
function of time (see Hopkins et al. 2006b for a discussion; 
the key point is that by Mb < -23, a more accurate predic- 
tion from our simulations of the HGLF may be up to a fac- 
tor ~ 2 lower), and as shown in Figure0the uncertainties in 
the initial MGLF (which propagate linearly into the HGLF) 
are large. For these reasons, the slight discrepancy in nor- 
malization between our prediction and the observations is not 
significant. Also, note that the appearance in Figure [Tsl that 
all bright mergers host optical quasars is an artifact of how 
steeply the Schechter function falls off - even at the brightest 
luminosities, there can be a factor r^2- 10 difference between 
the MGLF and Mb < -23 HGLF prediction. 

However, despite these uncertainties, the key qualitative 
point is clear - the quasar HGLF can be understood as a sub- 
set of the MGLF, with the appropriate though more sophisti- 
cated selection criteria applied as only those mergers during 
(or after) the "blowout" phase will be identified. The charac- 
teristic magnitude, relative normalization, and narrow width 
of the HGLF are explained in our modeling, demonstrating 
that this luminosity function is self-consistent with both the 
QLF and MGLF in our interpretation of quasar fueling. Ob- 
servations of e.g. the fraction of mergers with optical AGN 
(e.g.,|Sanchez et al. 2005), although traditionally more diffi- 
cult than observations of the fraction of AGN in mergers, also 
suggest a similar fraction of mergers hosting optical AGN to 
that we predict for the appropriate luminosity limits in Fig- 
ure^] Although we have only shown this at one redshift, 
where most of the observations have been made, the qualita- 



tive result is the same regardless of the redshift (or e.g. value 
of /gas or a) chosen. 

This does not imply that all quasar host galaxies will be 
visible as mergers or interacti ng systems, although the low- 
redshift sample of Bahcall et al. ( 1997) indeed identifies the 
vast majority as such. Generally, the "blowout" phase in our 
simulations follows the final coalescence of the black holes, 
meaning that tidal tails and other evidence of a recent merger 
should in principle be evident, but at higher redshift surface 
brightness dimming will make these features nearly impos- 
sible to observe. Moreover, once the blowout begins, star 
formation is terminated and the host galaxy rapidly begins to 
redden, and extended tidal features will fade in 1 Gyr 

The HGLF distribution we predict is consistent with ob- 
served properties of lower-luminosity AGN host galaxies. For 
example, Kauffmann et al. (2003a) study the host galaxies of 
~ 20,000 SDSS narrow line, relatively low-luminosity AGN 
at 7 < 0.3 and find that a large fraction of these objects corre- 
spond to what we expect from simulations for post-"blowout" 
objects with rapidly declining accretion rates in relaxing, 
rapidly reddening systems (see e.g. Hopkins et al. [2006a] for 
a description of the falloff in AGN luminosity as the merger 
relaxes). These objects reside in massive spheroids, with 
properties of "normal" ellipticals except for young stellar pop- 
ulations and evidence of starbursts in the past ~ 1 - 2 Gyr, with 
a sizable fraction ('^ 30%) of especially the brightest objects 
showing obvious evidence of interaction and/or recent merg- 
ers. Similar results are also found for e.g. the AGN host pop- 
ulation in GEMS ( Sanchez et al. 2004) and radio loud quasars 
("Sanchez & Gonzalez-Serrano 2003), although these may be 
preferentially at relatively low accretion rates (e.g.. .Ho..2002t 
iHopkins. Naravan. & Hernauistl2006l) . 

8. DISCUSSION 

We have used simulations of galaxy mergers which account 
for star formation, metal enrichment, radiative cooling, super- 
nova feedback and pressurization of a multi-phase interstellar 
medium, and black hole growth and feedback, to relate the 
distribution of observed quasar properties and its evolution 
with redshift to the distributions of merging galaxy proper- 
ties, including their luminosity and mass functions, character- 
istic gas fractions, and contribution to the star formation rate 
density of the Universe. Our simulations span a wide range 
of initial and final conditions, varying virial velocities, initial 
gas fractions, masses, orbital parameters, initial black hole 
masses, ISM gas equations of state, redshifts, and galaxy mass 
ratios. Our modeling allows us to self-consistently map be- 
tween the merging galaxy and quasar distributions in a phys- 
ically motivated manner without invoking tunable cosmolog- 
ical distributions and enables predictions of many properties 
of merging galaxies at different redshifts. 

We find that: 

• The joint scaling of dust obscuration and star formation 
with gas density yields quite flat lightcurves in the optical/UV, 
allowing us to provide a simple parameterization of the prob- 
ability for observing a particular merger at a certain luminos- 
ity. To enable future comparison with observations, we pro- 
vide fits from our simulations to infer the distribution of ob- 
served luminosities in different bands as a function of the stel- 
lar mass formed in mergers, total stellar mass, or final black 
hole mass of the merging systems (TableQ. Our simulations 
allow us to determine color and mass-to-light ratio distribu- 
tions of mergers, which we find to be in good agreement with 
those observed, and self-consistently map between e.g. ob- 
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served merger luminosity functions in the near IR, optical, 
and UV, star formation rate distributions, mass functions, and 
mass or luminosity-dependent merger rates. 

• We consider in detail observed merge r luminosity func- 
tion s in g-b and at z < 0.2 ("Xu et al.' 2004) and 280 nm at z = 
0.7 dWolf et al. 2005), and mass functions at z-- 0.4, 0.7, 1.2 
(iBundvetal. 2005a b). We use our simulations to map them 
to a quasar luminosity function at each redshift, and find that 
the predicted qua sar lumin osity functions agree with those ob- 
served by e.g. Ueda et al. (2003). Conversely, we invert this 
procedure, predicting the merger luminosity function from the 
observed quasar luminosity function, and again find agree- 
ment, although the predictions are significantly less well- 
constrained in this direction. Both distributions are self- 
consistent (each predicts the other) and can be used to predict 
merger rates and the rate of formation of quasars or spheroids 
as a function of that luminosity, mass and redshift. 

Although merger mass functions are not yet well de- 
termined above thes e redshifts, the demonstration in 
iHopkins et all j2006d) that the characteristic masses of ob- 
served mergers, the "quenching" or "transition" mass where 
elliptical galaxies begin to dominate the total galaxy popula- 
tion, and the characteristic mass of quasars (at the observed 
QLF break) all trace one another over the range < z < 3 
further suggests that this mapping is valid at high redshift. 
We extend our predictions of merger luminosity functions by 
using the observed QLF to predict the merger luminosity den- 
sity in K and B-bands from z = 0-6 with small factors ~ 2 
uncertainties. The prediction agrees wi t h existing observa- 
tions at z < 3 (Brinchmann et al. "1998'; Toledo et al. 19991 
IConselice et al. 2003, 2005; Xu et al. 2004; Bundv et alJ 
|2p05a b), but the present observational uncertainties are large. 
Future measurements of the mass or luminosity density in 
mergers will provide a critical test of this model and the asso- 
ciation between mergers and quasar triggering at all redshifts. 

•We compare this relation of merger and quasar distribu- 
tions to that inferred if quasars turn "on/off" in step-function 
fashion or follow exponential lightcurves, and find that the 
observations rule out such scenarios at > 99.9% confidence, 
even if we allow the quasar lifetime to vary freely to pro- 
duce the best fit. This can be determined entirely from the 
observed shapes of the QLF and MGLF, independent of se- 
lection efficiencies, gas fractions, or other cosmological prop- 
erties; it would require a > 4 - 6(T change in the observed 
shape to reverse this conclusion. Any comparison of merger 
rates or merger distributions and quasar activity must account 
for complex quasar lightcurves and the non-trivial nature of 
quasar lifetimes which vary as a function of both instanta- 
neous and peak luminosity. 

•We use the QLF to predict the distribution of star for- 
mation rates (SFRF) at various redshifts, and subsequently 
the SFR density in mergers as a function of redshift from 
z = - 6. Comparison of these predictions with measure- 
ments, which cur rently exist at z 5, 1 ( Brinchmann et al .ll998t 
iMenanteau et alJl2001[ l20()6t iBelTet al.u200 5). shows good 
agreement, and suggests that the QLF can provide a new, pow- 
erful, and independent probe of the SFR density in the Uni- 
verse driven by mergers. At redshifts z < 4, this allows us to 
predict the SFR density caused by mergers to a factor '--^ L5-2 
accuracy, with the primary systematic uncertainty being weak 
constraints on typical gas fractions of disks at low redshift 
z < L At z > 4, future improvements in constraints on the 
QLF break luminosity will enable similar predictions, with- 
out a large uncertainty in gas fractions, as gas-rich systems 



are expected at these high redshifts. Thus, quasar measure- 
ments can constrain the SFR density triggered by mergers at 
high redshifts and at the end of the reionization epoch, where 
attempts to measure morphological properties of galaxies and 
determine this quantity directly are not likely to be possible 
in the near future. This also provides a prediction of the total 
rate at which mass is moved from the blue (disk) sequence to 
the red (spheroid) sequence, independent of gas fraction, as 
analyzed in detail in Hopkins et al. (2006c). 

We find that the SFR density triggered in mergers is small at 
low redshift: ~ 1 % of the total SFR density at z 0. This frac- 
tion of the SFR density rises to a substantial, but not dominant 
~ 40% (upper limit 60%) at z ~ 2-3, and remains approxi- 
mately cons tant at higher redshift (rela tive to the total SFR 
density from Hopkins & Beacom ( 2006)). Given the low level 
of star formation in remnant ellipticals, the dominant contri- 
bution to the SFR density (especially at z < 2) must come 
from late-type spiral or irregular galaxies. 

• We predict the joint distribution of quasar and host 
galaxy luminosities in good agreement with observations 
tBahcall et al. 1997; McLure et al. 1999; Vanden Berk et afl 
120061) . The lack of a correlation between quasar and 
host galaxy luminosities (especia lly at low quasar lumi- 
nosity; Bahcall et al. 1997; McL ure et alJ 
'Jahnke & Wisotzki 2003; Heckman et al. 



2005; Vanden Berk et al. 2006) is a natural consequence of 
the nearly constant optical/U V merger light curves we discuss 
above while quasar luminosities change by orders of magni- 
tude during a merger. 

We predict the quasar host galaxy luminosity function 
(HGLF) from either the observed QLF or merging galaxy lu- 
minosity function (MGLF) at any given redshift. We find that 
the HGLF has a similar shape to the MGLF in a given band, 
but becomes increasingly peaked about M* as the minimum 
optical quasar luminosity or minimum ratio of optical quasar 
to optical host galaxy luminosity is increased. For a given 
minimum quasar optical luminosity, e.g. Mb < -23, we find 
good agreement between the observed HGLF ([Bahcall et'all 
1 19971 [Hamilton et al. 2002; Vanden Ber k et alJ20q6|) and our 
prediction from the observed MGLF of IXu et alJ i200 4^ or 
QLF of Ueda et al. (2003). The observed shape, normaliza- 
tion, and break/turnover luminosity in the quasar HGLF is a 
natural consequence of its being a subset of the MGLF, with 
e.g. a break at approximately the same luminosity (so long 
as the luminosity limit of the quasars in the sample is not so 
high as to exclude typical ^ hosts). This does not mean 
that quasar hosts will be obvious mergers, as the bright opti- 
cal phase of quasar activity is associated with the "blowout" 
phase following the final coalescence, in which difficult to 
observe, rapidly fading tidal features may represent the only 
morphological merger signature. 

We have demonstrated that the statistics and distributions 
of merging galaxy and quasar populations are self-consistent 
and can be used to predict one another in the context of 
the merger hypothesis. Coupled with the analysis of quasar 
properties in Hopkins et al. (2005a-d, 2006a,b), remnant red- 
sequence elliptical galaxies in Hopkins et al. (2006c), and 
the co-evolution of the "transition" mass and merger and 
quasar masses in Hopkins et al. (2006d), the modeling pre- 
sented here unifies the populations of multiple relevant merger 
stages. The number of observed mergers accounts for the 
bright quasar population as well as the observed buildup in 
the mass of the elliptical / red-sequence galaxy population, 
and vice versa. There is no room for a large fraction of gas- 
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rich mergers which do not produce a bright quasar phase and 
remnant, reddening elliptical, nor is there room for a large 
fraction of bright quasars or elliptical galaxies which are not 
formed in gas-rich mergers. 

This allows us to make a wide range of predictions and pro- 
vides a critical test of the hypothesis that starbursts, quasars, 
and elliptical galaxies are linked through the process of gas- 
rich mergers, a test which is complimentary and equally 
important to measurements of the individual photometric 
and kinematic properties and correlations of mergers and 
merger remnants (e.g., their profiles, metal enrichment, phase 
space densities, fundamental plane and Mm -a relations: 



Lake & Dressle Hll986t IPovon et alJll99l lO liva et all 
Shier & Fischei^ ll998t Uames et al.ll 19991 iGenzel et all 

■ ' ' ' '2006) ' 



122^ 



2001 



Rothberg & .Jose ph '2004, '2006) and comparison of these 
remnants to detailed simulations (e.g., Barnes & Her nquist 
[1992; Hernquist 1993a; Hernquist et al. 1993; Soringe l et all 
l2'005a; .Cox et al. .2005:. Robertson et al..2005a..2006b) . 

This does not mean, of course, that ellipticals never evolve 
via dry (gas-poor, spheroid-spheroid) mergers (see, e.g., 
iBell e t al. 2006a; van Dokkum 2005), as such mergers only 
mod ify a necessarily pre-exi sting population of ellipticals (see 
also iHopkins et alT i2006d ) for a calculation of the effects 
of such processes on observed elliptical distributions). This 
also does not imply that all AGN activity is associated with 
gas-rich mergers, as there can be a substantial contribution 
to low-luminosity AGN from activity triggered in quiescent 
spiral s (for a detailed comparison, see i Hopkins & Hernquist! 
120061) . as well as the large amount of quiescent, low-level ac- 
tivity in relaxed ellipticals which we predict (and account for) 
many dynamical times after the spheroid-forming merger as 
the quasar lightcurve decays (Hopkins et al. 2006b). 

Our modeling makes predictions that can be used to test 
our underlying theory. For example, mergers at high red- 
shifts should involve galaxies that are, on average, more gas- 
rich than local spirals (our preliminary comparison, for ex- 
ample, favors an exponentially d eclining f ^ ,^ ov er ^ a few 
Gyr). Preliminary evidence from iErb et al.l (12006 ) indicates 
that galaxies at redshift z~2 do, indeed, have large gas frac- 
tions /gas ^ 0.5, with some approaching /gas ^ 0.8-0.9, but 
future observatories such as ALMA should be able to measure 
this accurately. The connection between merger-driven star 
formation and quasar activity can also be used to test the cor- 
relation between our predicted SFR distributions and quasar 
peak luminosity. This should be possible with improved ob- 
servations of quasar hosts in relaxed, post-merger systems 
which can be used to determine their individual star formation 
histories (e.g.. Kauffmann et al. 2003a) . as well as measure- 
ments of star formation in obscured, lower-luminosity merger 
phases, possibly associated with IR-bright Type 2 quasar ac- 
tivity (e.g., C. Hao et al. 2005; Kim et al. 2006) and sti'ongly 
reddened optical quasars (e.g., Urrutia et al. 2005). 
Mergers of gas-rich spirals will imprint structure into the 



remnants that may be diffic ult to account for otherwise (see 
e.g., Robertson et al. 2006^ iCox et alJl2005l) . The starburst 
population left behind will characteristically modify the cen- 
tral light profiles and kinematics of merger remnants (e.g. Mi- 
hos & Hernquist 1994b; Hernquist & Barnes 1991), possibly 
explaining the central excesses of light seen in merging sys- 
tems (e.g. Rothberg & Joseph 2004, 2006), and the presence 
of kinematic subsystems in ellipticals. This can be tested in 
detail by comparing predictions for metallicity and color gra- 
dients and orbital distributions with observations. 

More subtle may be the shells, ripples, loops and other 
fine structures seen around many relaxed ellipticals (see, e.g. 
Schweizer 1998) that require a source of dynamically cold 
stellar material and are hence a natural consequence of ma- 
jor mergers involving disk galaxies (e.g. Hernquist & Spergel 
1992), but do not form in major mergers between hot stel- 
lar systems (e.g. Hernquist & Quinn 1988). A measurement 
of the rate of occurrence of fine structure in red galaxies with 
redshift would further constrain the importance of mer gers in- 
volving disk galaxies to the formation of ellipticals (iBell et alJ 
2006a; van Dokkum 2005). 

We use our simulations to develop a formalism to derive the 
relations between the various populations we have studied in 
a manner robust against different cosmological distributions 
which are poorly constrained and often tuned to reproduce ob- 
servations. However, this also means that we cannot constrain 
certain cosmological distributions. For example, we find that 
the distribution of observed luminosities of a merger is robust 
when expressed as a function of the final stellar mass or total 
stellar mass formed in the merger This enables us to map the 
MGLF to a merger mass function, but does not allow us to 
consider the relative contribution of mergers with e.g. differ- 
ent initial mass ratios, ISM equations of state, or star forma- 
tion histories. It will be interesting to see whether cosmologi- 
cal simulations yield the distributions of merger statistics that 
we have derived from quasar, merger, and elliptical galaxy 
observations. Coupling such cosmological descriptions with 
our detailed modeling of the complex star formation and black 
hole growth histories in mergers enables an a priori theoretical 
prediction of the wide array of phenomena we have demon- 
strated are linked through gas-rich galaxy mergers. 
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TABLE 1 

Best-Fit Observed Merger Luminosities 



Band 


Mm 


a 




M„ 


a 




Mg 
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X'l'^ 


cr 


280/40 


-20.7 ±0.1 


-2.0 ±0.1 


0.61 


-22.0 ± 0.2 


-2.0 ±0.3 


2.29 


-22.3 ±0.2 


-1.5 ±0.2 


2.76 


0.62 ± 0.06 


U 


-19.7 ±0.2 


-2.0 ±0.3 


0.76 


-20.9 ±0.2 


-2.3 ±0.3 


1.47 


-21.3 ±0.2 


-1.9 ±0.2 


1.59 


0.67 ±0.32 


B 


-20.2 ±0.2 


-2.0 ±0.2 


0.97 


-21.2±0.1 


-2.3 ±0.2 


0.83 


-21.7±0.1 


-1.7 ±0.2 


0.77 


0.49 ± 0.25 


V 


-20.7 ±0.1 


-2.0 ±0.2 


1.66 


-21.6±0.1 


-2.2 ±0.2 


0.74 


-22.1 ±0.1 


-1.7 ±0.2 


1.13 


0.34 ±0.28 


R 


-21.0±0.1 


-1.9±0.2 


2.40 


-21.8±0.1 


-2.2 ±0.2 


0.88 


-22.3 ±0.1 


-1.7 ±0.1 


1.61 


0.52 ±0.28 
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-21.2±0.1 


-1.9±0.1 


3.40 


-21.9±0.1 


-2.3 ±0.1 


0.97 


-22.5 ±0.1 


-1.7 ±0.1 


2.24 
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-22.2 ±0 A 


-2.0 ±0.1 


6.40 


-22.9 ±0.1 


-2.2 ±0.1 


1.51 
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3.69 


0.47 ±0.1 8 
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0.29 ±0.21 
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-20.7 ±0.2 
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Distribution of observed luminosities of mergers, as a function of new stellar mass M* ^ new formed during the merger. This assumes the probability of viewing a merger at some magnitude in the given band is a Gaussian with median 
(Mband) ^iiid dispersion cr in the band. 
^Same, as a function of total, final stellar mass M*^ ,oi after the merger. 
^Same, as a function of final black hole mass Mbh after the merger. 



